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features  of  the  surface  shear  stress. 
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ABSTRACT  OF  DISSERTATION 


A  TIME  DEPENDENT  FLOW  MODEL  FOR  THE  INNER  REGION 
OF  A  TURBULENT  BOUNDARY  LAYER 

Response  of  the  flow  variables  to  external  driving  forces  is 
non-linear  for  shear  flows.  For  the  turbulent  boundary  layer  case, 
surface  shear  stress  fluctuations  of  magnitude  as  great  as  the  mean 
value  are  observed.  For  flow  near  the  surface  Prandtl's  turbulent 
boundary  layer  approach  of  employing  averaged  Reynolds  equation  artd  a 
turbulence  closure  model  is  insufficient  to  account  for  surface  shear 
fluctuations.  A  model  which  incorporates  a  discrete  time  dependent 
solution  for  the  inner  region  of  the  turbulent  boundary  layer  is  pro¬ 
posed.  The  model  requires  stochastic  averaging  of  the  time  dependent 
solution  to  account  for  the  random  aspect  of  the  flow. 

The  physical  model  for  the  flow  near  the  surface  is  based  on  the 
bursting  cycle  observed  in  the  inner  region  of  a  turbulent  boundary 
layer.  Localized  pressure  gradients  created  in  ^ife  valleys  of  the 
large  scale  structures  of  the  outer  region  of  the  ^ low  are  assumed  to 
be  the  origin  of  the  bursting  process.  This  model  treats  the  sweep 
motion  as  an  impulsively  started  flow  over  a  flat  plate.  An  averaging 
technique  is  demonstrated  to  predict  the  important  features  of  the 
surface  shear  stress. 

In  order  to  confirm  the  time  dependent  model  assumptions, 
measurements  of  the  probability  distribution  and  cross-correlation  of 
the  longitudinal  turbulent  velocity  and  the  surface  shear  stress  were 
evaluated.  The  sweep-scale,  sweep-direction,  and  origin  of  the 
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instability  art'  <U‘Lt'rminc(l  from  isocorrolation  maps.  The  shape  of  the 
probability  density  distributions  of  the  velocity  near  the  surface  and 
the  surface  shear  stress  are  found  to  be  similar.  However,  the  velo¬ 
city  probability  distribution  changes  rapidly  with  increasing  distance 
from  the  surface. 

As  implied  by  the  time  dependent  model  for  the  surface  shear 
stress,  the  magnitude  of  the  large  surface  shear  stress  would  be  sub¬ 
stantially  changed  if  the  sweep  motion  could  be  modified.  A  series  of 
thin,  metal  plates  were  employed  to  block  the  instability  from  reaching 
the  surface.  Results  show  that  the  mean  value  of  surface  shear  and  the 
large  magnitude  fluctuations  of  surface  shear  stress  were  reduced 
significantly.  The  variation  in  surface  shear  was  found  to  be 
extremely  sensitive  to  slight  angle  of  attacks  of  the  plates. 
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CHAPTER  I 


INTRODUCTION 

Since  Prandtl  introduced  the  concept  of  a  boundary  layer  in  1904, 
the  viscous  effects  on  the  flow  adjacent  to  a  solid  boundary  have  re¬ 
ceived  a  great  deal  of  attention.  The  laminar  boundary  layer  problem 
has  been  solved  numerically  for  a  wide  range  of  flow  condition.  How¬ 
ever,  the  turbulent  boundary  layer  problem  is  still  far  from  being 
solved  due  to  its  complex  nature.  Attempts  which  paralleled  the  tech¬ 
niques  employed  in  solving  the  laminar  boundary  layer  problem  were  also 
applied  to  the  turbulent  cases.  The  eddy  viscosity  (or  mixing  length) 
concept  still  is  viewed  as  an  engineering  technique  to  evaluate  turbu¬ 
lent  shear  flows.  Conventionally,  a  model  which  divided  the  turbulent 
boundary  layer  into  an  inner  region  and  an  outer  region  was  hypothe¬ 
sized.  Considerable  amount  of  effort  has  been  made  to  obtain  better 
estimates  of  the  mixing  length  and  eddy  viscosity  for  the  two  distinct 
regions  of  a  turbulent  boundary  layer.  The  introduction  of  the  idea  of 
eddy  viscosity  substitutes  long-time  averaged  statistical  quantities 
for  the  time  dependent  properties  which  are  inherent  to  turbulent  flow. 
While  of  value  in  limited  engineering  applications,  the  early  models 
have  for  the  most  part  required  a  groat  deal  of  empirical  input. 

In  the  last  two  decades,  as  a  result  of  the  improvement  of 
experimental  techniques  and  the  advance  of  electronic  computer  tech¬ 
nology,  a  great  deal  of  information  about  the  turbulent  boundary  layer 
has  been  obtained.  The  new  information  has  led  to  a  better  understand¬ 
ing  about  the  structure  of  the  turbulent  boundary  layer.  Detailed  in¬ 
vestigation  of  flow  properties  in  the  near  wall  region,  and  the  outer 
region,  have  been  made.  Different  dominant  flow  features  exist  in  the 
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two  regiotts.  A  large  scale  motion  prevails  in  the  outer  region,  while 
the  bursting  phenomenon  is  the  most  important  feature  in  the  near  wall 
region.  However,  the  relationship  between  the  large  scale  motion  and 
the  bursting  phenomenon  is  not  well  understood.  It  is  also  necessary 
to  relate  the  features  to  the  important  properties,  such  as  surface 
shear  stress. 

In  the  present  study,  a  physical  model  which  describes  a  possible 
connection  between  the  large  scale  motion  and  the  bursting  phenomenon 
is  hypothesized.  Following  this  complex  model,  a  simplified  time  de¬ 
pendent  model  for  the  surface  shear  stress  under  a  turbulent  boundary 
layer  was  employed  to  illustrate  the  importance  of  time  dependent  solu¬ 
tion  of  a  turbulent  flow.  A  stochastic  averaging  technique  was  devel¬ 
oped  to  account  for  the  random  aspect  of  flow  in  predicting  the  surface 
shear  stress.  Predicted  surface  shear  was  compared  with  experimental 
results  by  assuming  a  Gaussian  or  a  modified  Rayleigh  probability  dis¬ 
tribution  of  the  sweep  motion.  Experimental  evidence  was  obtained  to 
support  the  time  dependent  model.  Experiments  were  performed  over  a 
nearly  zero  pressure  gradient  surface  in  a  small  wind  tunnel.  Time 
dependent  data  evaluated  include;  convective  velocity,  probability  and 
correlation  of  the  turbulent  velocity  and  surface  shear  stress.  A 
simple,  thin  plate  devices  was  used  to  modify  the  structure  of  the  flow 
near  the  surface,  which  in  turn  reduced  the  surface  shear  stress.  The 
device  was  based  on  the  implication  of  the  model  for  the  flow  in  the 
inner  region  and  its  relation  to  the  surface  shear  stress.  The  results 
appear  to  justify  the  assumed  model. 
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CHAPTER  II 


REVIEW  OF  THE  TURBULENT  BOUNDARY  LAYER  STRUCTURE 
2 . I  Early  Model 

Since  Prandtl  introduced  the  concept  of  a  boundary  layer  in  1904, 
the  viscous  effects  on  the  flow  adjacent  to  a  solid  boundary  have  re¬ 
ceived  a  great  deal  of  attention.  A  boundary  layer  could  be  either 
laminar  or  turbulent.  The  laminar  boundary  layer  problem  can  be  solved 
numerically,  however,  the  turbulent  boundary  layer  problem  is  still  far 
from  being  solved  due  to  its  complex  nature.  Pioneer  studies  made  by 
Prandtl,  von  Karman  and  others  contributed  considerably  to  the  early 
understanding  of  turbulent  boundary  layer  characteristics.  Convention¬ 
ally,  two  distinct  regions  in  a  turbulent  boundary  layer  were  hypoth- 
sized.  In  the  "wall",  or  "inner"  region,  the  viscous  effects  are 
important;  while  in  the  "outer"  region,  the  turbulent  transport  of 
momentum  is  dominant. 

Attempts  which  apply  the  same  technique  employed  in  solving  the 
laminar  boundary  layer  problem  were  also  used  to  evaluate  the  turbulent 
cases.  One  of  the  earliest  approaches  was  the  concept  of  eddy  visco¬ 
sity  introduced  by  Boussinesq.  Prandtl  constructed  the  mixing  length 
hypothesis  to  relate  the  turbulent  shear  term  to  the  mean  veocity 
gradient.  Followed  the  two-region  model  of  the  turbulent  boundary 
layer,  numerical  evaluation  of  the  mixing  length  assumed  that  it  in¬ 
creased  nearly  linear  with  distance  from  the  wall  in  the  wall  region, 
and  that  it  remained  nearly  constant  over  the  outer  region. 

Over  the  ensuing  years,  certain  refinements  of  the  model  were 
made,  mainly  regarding  the  evaluation  of  the  mixing  length  hypothesis. 
Van  Driest  (1956)  forsaw  the  fluctuating  nature  of  the  fluid  near  the 
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wall  and  suggested  a  damping  factor  for  the  eddy  viscosity  in  this 
area.  Townsend  (1956)  suggested  the  use  of  a  mixing  length  which  was 
corrected  in  accordance  with  the  intermittency  factor  in  the  outer 
region.  Other  refinements,  such  as  considering  the  entrainment  proper¬ 
ties  of  turbulent  boundary  layer,  Head  (1958),  have  led  to  some  im¬ 
provements  in  predicting  the  boundary  layer  properties. 

Parallel  to  the  above  model,  Clauser  (1956)  suggested  that  the 
outer  region  of  the  turbulent  boundary  layer  could  be  treated  as  a 
laminar  boundary  layer  having  a  thin  sublayer  of  a  different  fluid  with 
much  lower  viscosity  next  to  the  wall.  Measurements  of  flow  variables 
in  the  sublayer  region,  which  reflected  the  random  aspect  of  turbulent 
shear  flow,  made  before  1957  was  summarized  by  Corrsin  (1957).  These 
measurements  revealed  the  existence  of  large  magnitude  fluctuations  of 
surface  shear  stress,  surface  static  pressure  and  boundary  layer 
thickness . 

2.2  Detailed  Evaluation  of  the  Turbulent  Boundary  Layer  Structure 

In  the  last  decade  more  detailed  information  about  the  turbulent 
boundary  layer  structure  has  been  obtained.  This  new  information  has 
led  to  a  better  understanding  of  the  structure  of  turbulent  boundary 
layers;  however,  a  consistent,  workable  model  for  the  flow  has  proven 
elusive.  Recently  developed,  special,  photographic  techniques  have 
made  visualization  of  the  developing  turbulent  boundary  layer  possible. 
Figure  1  is  a  sketch  of  observations  and  Figure  2  shows  a  photograph  of 
a  turbulent  boundary  layer  developing  along  flat  plates  with  zero  pres¬ 
sure  reported  by  Falco  (1977)  and  Nagib  et  al.  (1979),  respectively. 

Different  aspects  of  the  flow  can  be  seen  in  the  two  figures. 
Figure  1  shows  the  overall  view  of  the  general  shape  of  the  large  scale 
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A  qualitative  picture  of  the  velocity  distribution  in  a  large 
scale  bulge  and  its  surrounding  fluid  was  obtained  by  Blackwelder  and 
Kovasznay  (1972),  using  conditional  sampling  techniques  to  evaluate  the 
hot  wire  signals.  Figure  3  shows  their  results,  wherein  U  =  0.93  U 
represents  the  mean  velocity  of  the  turbulent  bulge.  The  rotational 
nature  of  the  large  scale  motion  in  the  bulge  was  demonstrated  by 
Blackwelder  and  Kovasznay  (1972).  The  local  velocities  at  the  front 
and  back  of  the  interfaces,  measured  by  Kaplan  and  Laufer  (1969)  using 
a  ten-hot-wire  rake  across  the  boundary  layer,  are  shown  in  Figure  4. 
It  was  found  that  the  downstream  side  of  the  interface  moves  faster 
than  the  upstream  side.  Kovasznay  et  al.  (1970)  reported  zone  averages 
of  the  fluctuating  and  streamwise  mean  velocity  components,  as  shown  in 
Figures  5  and  6.  Within  the  non-turbulent  zone  the  fluid  moves  faster 
than  it  does  in  the  turbulent  zone.  The  turbulent  intensity  in  the 
non-turbulent  region  is  lower  than  in  the  turbulent  zone,  though  it  was 
not  zero. 


Figure  3.  Velocity  distribution  in  the  outer  region  of  the  boundary 
layer  obtained  by  Blackwelder  and  Kovasznay  (1972). 
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Figure  4. 


Velocity  distribution  at  the  turbulent  bulge  interfaces 
obtained  by  Kaplan  and  Laufer  (1969). 
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Figure  5.  Zone  averages  of  the 
streamwise  velocity  component 
reported  by  Kovasznay  et  al. 
(1970). 


Figure  6.  Zone  averages  of  the 
intensity  of  the  streamwise  ve¬ 
locity  fluctuations  reported  by 
Kovasznay  et  al.  (1970). 


From  the  illustrations  jt  Figures  3,  4  and  5,  it  appears  that  in 
the  valleys  of  the  consecutive  large  eddy  structures  there  is  a 


localized  pressure  gradient,  which  tends  to  push  the  fluid  inward 
toward  the  wall,  and  also  accelerates  the  fluid  on  the  upstream  side  of 
the  turbulent  bulge.  The  local  streamline  curvature  near  the  wall, 
which  would  imply  a  local  pressure  variation  is  indicated  in  Figure  2. 
The  non-turbulent  fluid  is  seen  to  be  thrust  almost  to  the  surface. 

Although  the  photographs  made  by  Falco  (1977)  and  Nagib  et  al. 
(1979)  show  the  large  scale  coherent  structure  in  the  outer  region  of 
the  turbulent  boundary  layer,  the  structure  buried  in  the  confined  wall 
region  is  not  discernible.  Special  equipment  and  techniques  are  being 
used  to  explore  the  turbulent  structure  in  the  wall  region.  Kline 
et  al.  (1967)  conducted  visual  studies  by  using  a  hydrogen  bubble  tech¬ 
nique.  Corino  and  Brodkey  (1969)  observed  motions  of  suspended  colloi¬ 
dal  particles  in  the  wall  region  of  a  circular  pipe  flow  by  using  a 
high-speed  camera  moving  with  the  flow.  Combined  visualization  and 
hot-wire  anemometer  techniques  were  also  employed  by  researchers  such 
as  Kim  et  al.  (1971),  Falco  (1977,  1980)  and  others. 

Visual  studies  of  Kline  et  al.  (1967),  Corino  and  Brodkey  (1969) 
and  Kim  et  al.  (1971)  revealed  the  existence  of  a  somewhat  well- 
organized  but  spatially  and  temporally  dependent  motions  within  the 
wall  region.  A  large  scale,  streaky  structure  and  an  intermittently 
occurring,  violent  bursting  process  were  observed  in  the  wall  region. 
A  cycle  of  the  bursting  phenomenon  was  described  with  the  help  of 
sketches,  Figure  7,  by  Corino  and  Brodkey  (1969).  Each  sketch  of 
Figure  7  shows  an  important  step  in  the  motion  involved  in  a  bursting 
cycle. 

First,  near  the  wall  low-speed  streaks  form  and  gradually  grow  to 
a  vertical  dimension  of  y*  ^  10,  as  shown  in  Figure  7a).  When  the 


streak  reaches  y*  ~  10  it  starts  to  oscillate  and  continues  to  grow. 
The  oscillation  amplifies,  as  it  continue  to  rise,  until  it  becomes 
unstable  and  suddenly  breaks  into  turbulent  motion  at  a  height  of  10  < 
y*  <  30.  Simultaneously,  a  portion  of  the  low  speed  fluid  of  the 
streak  is  ejected  into  the  core  region,  which  is  termed  the  "ejection" 
motion,  and  shown  in  Figure  7b).  Following  the  ejection,  a  mass  of 
fluid,  with  dimensions  larger  than  the  ejection  scale  and  with  velocity 
greater  than  the  local  mean  velocity,  flows  parallel  to  or  at  a  slight 
angle  toward  the  wall.  The  inward  flow  of  the  high-speed  fluid  is 
called  the  "sweep"  and  i.s  shown  in  Figure  7c).  Kim  et  al.  (1971) 


Figure  7.  Steps  in  the  flow  cycle  near  the  wall.  Reported  by 
Corino  and  Brodkey  (1969). 
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provided  photographs  of  streak-lines  using  hydrogen  bubbles  for  each 
stage  of  the  bursting  phenomenon.  Readers  are  referred  to  Kim's  report 
for  more  details. 

Associated  with  the  sequence  of  motions  in  bursting,  occasionally 
part  of  the  low  speed  ejected  fluid  is  deflected  back  toward  the  wall 
and  at  the  same  time  incoming,  accelerated  fluid  is  reflected  outward 
away  from  the  wall.  These  motions  are  called  "inward  interaction"  and 
"outward  interaction"  respectively.  Combinations  of  streamwise  and 
vertical  velocity  fluctuations,  u'  and  v' ,  as  shown  in  Table  1  are  used 
to  designate  these  motions. 


Table  1.  Signs  of  u' ,  v'  and  u'v*  related  to  motions  in  a 
bursting  cycle. 


Type  of  Motion 

Sign  of  u’ 

Sign  of  v’ 

Sign  of  u'v' 

Ejection 

- 

+ 

- 

Sweep 

+ 

- 

- 

Inward  Interaction 

- 

- 

•f 

Outward  Interaction 

+ 

+ 

+ 

Flow  properties  of  each  motion,  especially  the  Reynolds  shear 
stress  -  p  u'v',  measured  in  the  near  wall  region  by  Willmarth  and  Lu 
(1971),  and  Brodkey  et  al.  (1974),  indicate  the  importance  of  bursting 
in  the  generation  of  turbulent  energy.  Results  of  Willmarth  and  Lu 
(1971)  showed  that  fractional  contribution  to  -  p  u'v'  at  y*  =  30  are 
80%  and  43%  due  to  the  ejection  and  sweep  respectively. 

Although  the  bursting  phenomenon  is  spatially  and  temporally 
unsteady,  agreement  on  instantaneous  traverse  streak  spacing  and  mean 
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bursting  period  from  different  measurements  was  obtained.  Results 

obtained  in  water  channels  (1800  <  <  2500)  by  Kline  et  al.  (1967), 

Bakewell  and  Lumley  (1967),  Kim  et  al.  (1971)  and  Gupta  et  al.  (1971) 

showed  a  nondimensional  streak  spacing  of  z*  ~  100.  Rao  et  al.  (1971) 

summarized  the  results  of  Kim  et  al.  (1971),  Runstadter  et  al.  (1963), 

Laufer  and  Badri  Narayanan  (1971),  and  shown  that  the  mean  burst  time 

period,  T,  scales  with  the  free  stream  velocity,  U^,  and  boundary  layer 

thickness,  6;  and  TU^/6  ~  5  in  the  ran.-e  500  <  R^  <  9000.  Results 

reported  by  Falco  (1980)  shown  that  fu  /6  could  range  from  4  to  10  for 

00  ” 

approximately  the  same  Rg  range. 

As  demonstrated  above,  the  details  of  the  bursting  process  in  the 
wall  region  is  well  documented.  The  consequent  coherent  structure  in 
this  region,  on  the  other  hand,  has  not  received  equivalent  attention. 
Bakewell  and  Lumley  (1967)  and  Lee,  Eckelman  and  Hanratty  (1974)  were 
among  those  who  proposed  streamwise,  concentrated  vortices  in  the  wall 
region.  Recent  measurements  of  the  cross-correlation  between  the  sur¬ 
face  shear  stress  and  the  velocity  field  in  the  inner  region,  made  by 
Kreplin  and  Eckelman  (1979),  have  confirmed  the  existence  of  streamwise 
vorticity  in  the  sublayer.  The  measurements  suggested  that  the  vorti- 
city  exists  as  a  counterrotating  vortex  pair.  In  this  regard,  Falco 
(1980)  presented  results  which  he  obtained  by  using  visualization 
techniques  and  a  hot-wire  anemometer  simultaneously,  which  confirm  the 
presence  of  the  vortex  motion. 

Falco  applied  oil-fog  contaminant  through  a  slit  in  the  wall  under 
a  turbulent  boundary  layer  and  took  pictures  of  the  patterns  of  the 
oil-fog,  which  is  carried  by  the  coherent  motion.  Figure  8  is  a  photo¬ 
graph  reported  by  Falco  (1980).  The  pocket  structure  where  the  smoke 
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was  washed-out,  illustrates  the  existence  of  slreamwise  vortices  in  the 
near  wall  region.  Falco  djvided  the  evolution  of  a  pocket  into  five 
stages.  The  flow  properties  sucli  as  u,  du/dy  and  uv  for  each  stage 
were  evaluated.  Connection  of  each  stage  to  the  bursting  process  was 
discussed  in  great  detail.  Falco  concluded  that  there  exists  a  turbu¬ 
lence  generation  mechanism  near  the  will.  This  mechanism  includes, 
stretching  of  the  vorticity  in  tin  sweep  motion,  the  generation  of 
vorticity  near  the  wall  by  tin-  ;;lagn.ilion  point  flow  that  the  sweep 
creates  and  mutual  interaction  cl  the  vorliies  f mined  as  a  result  of 
these  processes  leading  to  motion  away  I  rein  tlie  wall. 


POCKETS 


Figure  8.  Sketch  iif  the  pocket  structure  near  tlie  wall  photographed 
by  Falco  (1930)  jt  =  7 '38. 

Several  other  flow  characteristics,  which  are  directly  related  to 
turbulent  boundary  layer  structure,  have  also  been  investigated.  The 
surface  pressure  variatuns  under  turbuient  boundary  layers  have  been 
studied  in  detail.  Ro^-.oa  rchc  r  s ,  such  as  Corcos  (196o),  Blake  (1970), 
Willmarth  and  Rons  (1965)  ind  others  measured  spai e-lime  correlations 
between  fluctuations  of  surface  pressure  and  turbulent  velocities.  The 
convective  nature  of  the  wall  pressure  was  studied  rather  extensively. 
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Favre  et  al.  (1957),  Wills  (1964,  1967)  and  others  measured  space-time 
correlations  of  turbulent  velocities  at  two  locations  with  the  stream- 
wise  separation  across  the  boundary  layer.  Their  results  showed  the 
frequency  dependent  nature  of  the  convective  velocity.  Also,  in  the 
outer  region  the  convective  velocity  was  smaller  than  the  mean  flow 
velocity  and  larger  in  the  wall  region.  Cliff  and  Sandborn  (1973) 
proposed  a  physical  model  for  the  convective  velocity  in  a  turbulent 
boundary  layer,  which  also  predicted  this  trend.  The  model  of  Cliff 
and  Sandborn  (1973)  postulated  that  packets  of  turbulent  fluid  are 
generated  in  a  production  zone  near  the  viscous  sublayer.  These 
packets  were  found  to  be  discernible  from  the  mean  motion  and  may  move 
either  outward  from  or  inward  to  the  wall.  The  magnitude  of  the  con¬ 
vective  velocity  in  the  production  zone  was  found  to  be  equal  to  the 
local  mean  velocity. 

Large  magnitude  fluctuations  in  surface  shear  stress  were  first 
mr'asured  by  Mitchell  and  Hanratty  (1966)  using  an  electrochemical 
technique.  Blirico  and  Simons  (1974),  and  Sandborn  (1979)  also  reported 
large  fluctuations  of  surface  shear  stress  both  in  water  channels  and 
in  wind  tunnels.  These  results  demonstrated  that  the  surface  shear 
stress  fluctuations  in  turbulent  boundary  layer  flows  are  highly  time 
dependent  and  of  large  magnitude.  The  large  fluctuation  characteris¬ 
tics  of  the  surface  shear  stress  have  been  used  as  a  means  of  studying 
the  turbulent  layer  structure.  Brown  and  Thomas  (1977)  reported  a 
limited  number  of  correlation  measurements  between  the  surface  shear 
stress  and  the  turbulent  flow  field.  The  surface  shear  fluctuations 
were  employed  by  Kreplin  and  Eckelmann  (1979)  for  studies  regarding  the 
tr.insverse  spacing  of  the  low-speed  streaks  and  streamwise  vortices. 


Based  on  the  above  observations,  the  time  dependent  nature  of  the 
turbulent  boundary  layer  becomes  evident.  Since  the  introduction  of 
the  idea  of  eddy  viscosity,  boundary  layer  models  have  tended  to  sub¬ 
stitute  long-time  averaged  statistical  quantities  for  time  dependent 
properties.  By  doing  so,  the  real  mechanisms  inherent  to  a  turbulent 
boundary  flow  often  have  been  overlooked  and  the  progress  toward  their 
understanding  limited.  The  response  of  a  boundary  layer  type  shear 
flow  to  a  change  in  velocity  or  pressure  is  nonlinear.  The  statistical 
averaged  approach  to  the  problem  masks  the  nonlinear  effects,  such  as 
the  observed  large  variations  in  surface  shear  stress.  With  the  obser¬ 
vation  of  coherent,  repeatable,  structures  within  the  boundary  layer, 
time  dependent  models  for  the  flow  should  be  possible.  Solution  of  the 
time  dependent  problem  first,  followed  by  the  use  of  a  statistical 
averaging  technique  should  help  retain  the  nonlinear  aspects  of  the 
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SURFACE  SHEAR  STRESS  FLUCTUATIONS 

To  predict,  the  flow  variables  in  the  inner  region  of  a  turbulent 
boundary  layer,  where  large  magniture  fluctuations  have  been  observed, 
the  non-linear  interaction  between  the  flow  parameters  need  be  con¬ 
sidered.  Based  on  the  reported  experimental  results,  a  physical  model 
for  the  flow  near  the  surface  is  hypothesized.  A  time  dependent  solu¬ 
tion  for  a  simplified  versiofi  of  this  model  was  pursued.  Utilizing 
this  solution  and  assuming  a  velocity  distribution  for  the  sweep  motion 
a  stochastic  averaging  technique  was  developed  to  predict  the  probabil¬ 
ity  density  distribution  of  surface  shear  stress.  The  purpose  of  this 
study  is  to  illustrate  the  importance  of  the  philosophy  required  to 
deal  with  large  magnitude  fluctuations  in  a  turbulent  shear  flow.  The 
philosophy  requires  that  the  time  dependent  equations  of  motion  be 
solved  first  and  then  a  statistical  average  be  employed  to  account  for 
the  random  aspects  of  the  turbulent  flow. 

3. 1  Physical  Model 

As  discussed  in  the  last  chapter,  visual  studies,  made  by  Kline  et 
al.  (1967),  Corino  and  Brodkey  (1969),  Kim  et  al  .  (1971)  and  others, 
revealed  that  in  the  sweep  motion  associated  with  the  bursting  phenom¬ 
enon,  a  stream  of  high  speed  fluid  from  the  outer  region  was  pushed 
into  the  wall  region.  This  high-speed  flow  will  create  an  impulsive 
effect  on  the  flow  in  the  proximity  of  the  wall.  It  appears  reasonable 
that  the  sweep  motion  is  responsible  for  the  appearance  of  large  magni¬ 
tude  surface  shear  stress  fluctuations,  as  observed  by  Mitchell  and 
Hanratty  (1966),  Blinco  and  Simons  (1974),  and  Sandborn  (1979). 
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As  first  argued  by  Corrsin  (1957),  tt'is  large  magnitude 
fluctuation  of  surface  shear  stress  would  not  be  seen  if  the  time 
averaged  boundary  layer  equations  were  solved.  Due  to  the  nonlinear 
relation  between  the  external  drivjng  force  and  tlie  flow  response  in  a 
turbulent  boundary  layer,  time  averaging  of  the  driving  force  removes 
the  extremes  of  the  flow  response.  Thus,  it  is  evident  that  a  Linif 
dependent  approach  is  needed  to  account  for  this  rion-linear  response 
characteristic  of  viscous  flows. 

Relating  the  surface  shear  stress  to  the  sweep  motion,  together 
with  the  observations  of  Rao  et  al.  (1971)  and  Falco  (1980)  that  the 
mean  bursting  period  scales  with  the  outer  region  parameters,  a  physi¬ 
cal  model  for  the  surface  shear  stress  in  a  turbulent  boundary  layer 
can  be  hypothesized. 

Near  the  surface  in  the  valleys  between  the  large  eddies 
associated  with  the  outer  region,  there  exists  strong,  localized, 
vertical  pressure  gradients.  At  the  beginning  of  a  sweep  sequence, 
there  is  a  high  shear  between  the  retarded  low-speed  streak  and  the 
outer  region  turbulent  flow.  When  the  suddenly  accelerated  fluid 
associated  with  the  sweep  motion  intrudes  onto  the  valley,  the  shear 
increases  and  a  large  scale  vortex  motion  appears.  The  flow  then 
breaks  into  turbulent  motion  and  packets  of  traverse  vorticity  are 
formed.  The  vorticity  is  carried  into  the  large  eddy  motion  by  the 
surrounding  velocity  field.  With  the  help  of  the  pressure  gradients  in 
the  valleys  between  the  large  eddies,  a  stream  of  outer  flow  migrates 
inward  momentarily  to  meet  the  requirement  of  continuity.  The  downwar<l 
high-speed  flow  forms  streamwise  vortices  with  the  surrounding  low 
velocity  field  and  moves  downstream.  When  these  .streamwise  vortices 
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impulsively  arrive  in  the  proximity  of  the  wall,  very  high  magnitudes 
of  surface  shear  stress  are  created.  The  viscous  effect  retards  the 
high-speed  fluid  near  the  wall  in  a  Rayleigh  like  suddenly  accelerated 
plate  fashion.  The  retardation  of  flow  results  in  the  formation  of  the 
low-speed  streaks. 

The  visualization  results  ol  Magib  et  al.  (1979),  as  shown  in 
Figure  2,  together  with  the  flow  pattern  observed  by  BlackweTder  and 
Kovasznay  (1972),  Figure  3,  provided  evidence  for  the  existence  of  a 
localized  pressure  gradient.  As  for  the  sequence  of  the  ejection  and 
the  sweep  motions,  visual  studies  of  Kline  et  al  .  (1967),  and  Corino 
and  Brodkey  (1969)  indicated  that  the  sweep  motion  occurs  after  the 
ejection  motion  but  an  ejection  does  not  necessarily  cause  a  sweep 
motion.  Falco  (1977)  also  observed  this  phenomenon.  Falco  reported 
hot-wire  measurements  of  u’ ,  v’  and  u'v*  at  y'"  =  67  and  two  types  of 
large  scale  motion  measured  at  this  location,  as  shown  in  Figures  9a) 
and  9b).  He  explained  that  sometimes  the  motion  shown  in  Figure  9b) 
follows  that  in  Figure  9a),  but  a  regular  pattern  is  not  apparent. 
This  also  explained  the  intermittent  nature  of  the  bursting  process. 
Falco  also  implied  that  the  sweep  motion  follows  the  ejection  motion 
but  not  vise  versa. 

The  penetration  of  the  high-speed  outer  flow  into  the  sublayer  can 
be  demonstrated  by  the  following  measurements.  Chen  and  Blackv.elder 
(1978)  used  temperature  as  a  passive  contaminant  in  the  inner  flow.  By 
heating  the  entire  wall  of  a  wind  tunnel  test  section  to  approximately 
12°C  above  the  free-stream  temperature,  a  sharp  internal  temperature 
front,  characterized  by  a  rapid  decrease  in  temperature,  was  found  to 
extend  throughout  the  boundary  layer.  Their  results  are  illustrated  in 
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Figure  9.  Ensemble  averaged  u' ,  v'  and  u'v'  signals  for  large  scale 
motion  measured  at  y*  =  67  by  Falco  11977). 

a)  Large  scale  motion  type  1  has  a  negative  zone  averaged 
streamwise  velocity  perturbation. 

b)  Large  scale  motion  type  2  has  a  positive  zone  averaged 
streamwise  velocity  perturbation. 

Figure  10,  where  the  arrows  indicate  a  temperature  front.  Since  the 
internal  temperature  front  is  an  indic.ation  of  the  upstream  side  of  the 
large  eddy  structure  (back  of  large  eddy  as  used  in  Figures  3  and  4), 
the  data  indicated  that  the  sharp  acceleration  associated  with  the 
bursting  phenomenon  does  indeed  penetrate  down  to  the  sublayer. 
Ecklemann  (1974)  placed  hot-film  probes  right  above  a  surface  film  and 
recorded  signals  representing  the  streamwise  velocity  and  the  surface 
shear  stress  simultaneously.  Figure  11  shows  a  set  of  these  results 
for  y*  =  1.0,  1.9,  2.9,  4.8  and  8.6.  An  almost  one-to-one  correspon¬ 
dence  between  the  streamwise  velocity  in  the  proximity  of  the  boundary 
and  the  surface  shear  stress  is  noted.  This  evidence  is  used  to  jus¬ 
tify  the  belief  that  the  surface  shear  stress  is  directly  related  to 


the  sweep  motion. 
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Figure  10.  Simultaneous  temperature  .signal.s  from  the  ten-wire  rake 
in  the  turbulent  region  reported  by  Chen  and  Blackwelder 
(1978).  One  particular  temperature  front  is  denoted  by 
the  arrows. 
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Figure  11.  Simu I t.)neous  time  records  of  instantaneous  (3u/9y)o  and 
u  fluctuations  .at  v.irious  y*  position  reported  by 
F.ckelmann  (1974). 
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3.2  Simplified  Time  Dependent  Model 

Direct  conversion  of  the  model  of  a  sweep-burst  motion  into  an 
analytical  solution  for  the  turbulent  boundary  layer  could  be  extremely 
complex.  In  order  to  develop  an  analytical  approach  it  was  necessary 
to  evolve  a  simplified  model.  The  basic  assumptions  for  the  simplified 
model  are  as  follows: 

(a)  The  time  dependent  boundary  l.iyer  approximation  lor  the  Navier- 
Stokes  equations  for  a  two-dimensional  flow  over  a  flat  plate  is 
assumed  to  apply  in  the  sublayer  region. 

(b)  The  time  dependent  solutio/i  starts  with  the  beginning  of  the  sweej) 
motion  portion  of  the  overall  sweep-burst  cycle. 

(c)  The  localized  vertical  pressure  gradient  across  the  sublayer, 
which  must  in  part  initiate  the  sweep  motion,  is  of  a  very  short 
time  duration.  The  vertical  pressure  gradient  is  assumed  to 
disappear  instantaneously  after  the  sweep  motion  starts. 

(d)  The  initial,  equivalent,  free-stream  velocity  is  that  of  the 
non-turbulent  fluid  which  penetrates  into  the  sublayer  region  .at 
the  beginning  of  the  sweep  motion.  Its  magnitude  is  expected  to 
correspond  to  the  values  obtained  by  extrapolation  of  the  measure¬ 
ments  (non-turbulent  velocity)  to  y*  ~  30  of  Kovasznay  et  al . 
(1970),  as  shown  in  figure  5.  The  equivalent  free-stream  velocity 
should  also  be  related  to  the  convective  velocity  in  the  sublayer 
region . 

For  the  case  of  zero  pressure  gradient  in  the  longitudinal 
direction,  the  two  dimensional,  time  dependent  equation  of  motion  in 
the  x-di reel  ion  is: 


9u 

9t 
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Du  ^  Du 
9x  ''  Dy 


D^i 


9x 


(1) 


Using  the  results  of  convective  velocity  and  cross-correlation 
measurements,  which  are  discussed  in  Chapter  V,  the  above  equation  can 
be  simplified  by  estimating  the  order  of  magnitude  for  each  term.  For 
the  case  of  =  10.20  m/sec,  the  sweep  motions  originate  approximately 
at  0.7  cm  above  the  boundary  and  7.0  cm  upstream  of  the  surface  shear 
sensor.  The  averaged  period  of  occurrence  of  sweep  motion  was  esti¬ 
mated  as  28  milliseconds,  which  was  the  time  delay  for  maximum  correla¬ 
tion  between  a  hot-wire  signal  placed  at  the  above  coordinate,  (7.0, 
0.7),  and  the  surface  shear  sensor.  The  characteristic  velocity  asso¬ 
ciated  with  the  sweep  motion  was  estimated  as  the  convective  velocity 
measured  at  the  height  where  it  has  the  same  magnitude  as  the  local 
mean  velocity.  u  S  5.0  m/sec  was  selected  based  on  the  measurements. 
By  assuming  the  longitudinal  velocity  is  an  order  of  magnitude  larger 
than  the  vertical  velocity  within  the  flow  domain  considered,  the  order 
of  magnitude  for  each  term  can  be  estimated.  The  presence  of  a  verti¬ 
cal  pressure  gradient  and  the  ejection  motion  may  produce  locally  large 
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It  is  found  that  the  term  v  9  u/9x  is  much  smaller  than  the  other 
terms,  thus  it  can  be  neglected  similar  to  the  results  for  the  time 
averaged  equations.  The  resulting  Navier-Stokes  equation  for  ‘he  time 


dependent  flow  in  the  x-direction  is 
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The  truncation  of  the  time  dependent  solution  of  Equation  (2)  can  be 
estimated  as  equivalent  to  T,  which  is  the  most  probable  period  between 
two  consecutive  bursting  cycles. 

In  the  present  study,  the  model  was  developed  to  predict  the  time 
history  of  the  surface  shear  stress  during  the  bursting  cycle.  The 

flow  domain  considered  is  the  inner  region  of  the  boundary  layer.  A 
variation  of  the  effective  freestream  velocity  might  be  expected  over  a 
bursting  cycle  duration,  which  also  perturbs  the  flow;  however,  this 

effect  should  be  small.  The  response  of  the  surfact  shear  stress  to  a 
variation  of  freestream  velocity  has  been  investigated  by  Watson  (1956) 
for  the  case  of  a  laminar  boundary  layer  with  surface  suction.  The 
surface  shear  stress  response  was  demonstrated  to  be  much  smaller  than 
that  caused  by  an  impulsively  started  motion.  Thus,  it  was  assumed 
that  small  variations  in  the  effective  freestream  velocity  would  be  of 
secondary  importance. 

There  are  certain  limitations  imposed  on  the  model.  Since  the 

ejection  and  the  sweep  motions  are  large  scale  vertical  motions,  re¬ 

finement  of  the  model  could  require  inclusion  of  the  equation  of  motion 
in  the  y-direction  also.  Assumption  (c)  can  be  taken  to  imply  that  a 
zero  thickness  boundary  layer  exists  instantaneously  as  the  sweep 
motion  starts.  Thus,  an  infinitely  large  surface  shear  stress  would  be 
required  by  the  model.  Obviously,  it  will  be  necessary  to  limit  the 
surface  shear  to  finite  values. 

As  the  bursting  process  is  very  spatially  and  to  a  lesser  extent 


temporally  dependent,  the  strength  of  each  sweep  motion  past  a  given 


V,  \ 


point  will  vary.  A  complete  solution  of  the  mean  flow  properties  will 
require  a  stochastic  averaging  of  the  individual  sweep-burst  cycles. 

Consider  the  surface  shear  stress  induced  by  the  "sweep"  motion  in 
the  near  wall  region  of  a  turbulent  boundary  layer.  Assume  that  as 
soon  as  the  "jet-like"  stream  carried  by  the  sweep  reaches  the  wall,  a 
flow  field  is  established  instantaneously.  Then,  for  the  case  of  zero 
pressure  gradient,  the  governing  equation.  Equation  (2)  and  the  con¬ 
tinuity  equation  (9u/9x  +  9v/9y  =  0)  are  solved.  The  boundary 
conditions  are; 


for  t  >  0 ,  y  ^  u  =  -■  constant  (3) 

for  t  >  0 ,  y  =  0 ,  u  =  0 ,  v  =  0  (4) 


which  is  the  unsteady  Blasius  problem. 

By  introducing  the  following  nondimensional  terms, 
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u  =  u/U  ;  X  =  x/L;  t  =  -j— 
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Equations  (1)  together  with  the  continuity  equation  and  the  boundary 
conditions  become: 
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t°  > 

0,  y°  ^  ®, 

u  =  1 

(8) 

t°  > 

0,  y"  =  0, 

0  On 

U  =  V  =  0 

(9) 

The  above  equations  are  further  transformed  by  using 
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=  —  and  n  =  y°/^x° 


(10) 


as  the  independent  variables  and 


^  o  j  o  o  /  i- 

f  =  u  and  g  =  V  y  /r)  -  fr| 


The  transformed  equations  become; 
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and 
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If  t*  =  1/(1)  =  t'^/x°  is  introduced,  Equations  (12)  and  (13)  become: 


-2^’’'  Ip  *  In  "  '  -  “ 


(2  .  2t*f)  2£_  ,  ,  y  .  li  ,  0 

an 


The  Transformed  boundary  conditions  are: 


(14) 


(13) 


g  =  0  at  n  =  0 


and 


f  =  1  at  n 


A  numerical  solution  for  Equations  (14)  and  (15)  with  the  same  boundary 
conditions  was  reported  by  Watkins  (1*175)  for  unsteady  heat  transfer  in 
impulsive  Falkner-Skan  flows.  An  approximate  soluLioti  for  the  same 
problem  was  later  given  by  Gotti fredi  and  Quiroga  (1978).  In  both 
studies,  calculations  of  the  non-dimensional  surface  shear  stress  for 
several  flow  conditions  were  also  presented.  A  comparison  of  the 
results  obtained  from  the  two  solutions  were  shown  to  be  in  good 
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agreement.  Gottifredi  and  Quiroga  obtained  an  approximate  form  of  the 
non-dimensional  surface  shear  stress  for  the  general  Falkner-Skan 


flows.  Unfortunately,  an  error  was  found  in  this  approximate  form. 
Their  techniques  were  applied  to  the  present  problem  and  the  correct 
form  was  employed. 

By  defining  a  new  variable  ^  as 

4  =  n/V2t*  (16) 

and  assuming  a  series  solution  for  f  and  g  of  the  following  forms: 

f  =  Fq(4)  +  Fj(4)t*  +  ....  (17) 

g  =  Gj(4)t*^/^  +  G^(C)t*^^^  +  ...  (18) 

Equations  (14)  and  (15)  become  two  approximate  formulas  which  contain 
terms  involving  Fq(4),  Fj(4),  G^(4)  and  62(4) •  By  collecting  terms  of 
like  power  of  t*  the  following  ordinary  differential  equations  were 
obtained . 


24Fo’ 
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2eFj' 

*  S"  -  =  "Vo’  SS’ 

(20) 

S’ 

V2"4Fq'  t  V2Fq  =  0 

(21) 

Detailed  derivations  which  lead  to  the  above  equations  are  given  in 
Appendix  A.  Boundary  conditions  for  Equations  (17)  and  (18)  become: 

=  F,  =  G,  =  0  .t  t  =  0 
Fp  =  1  and  F^  =  0  at  4  “ 

To  solve  Equations  (19),  (20)  and  (21),  an  approximate  technique  was 

used,  which  was  attributed  to  Rosenzweig  (1959)  (and  was  applied  by 
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Cess  (1961)  to  the  problem  of  heat  transfer  due  to  a  nonsteady  surface 
temperature).  This  technique  requires  an  estimate  of  df/d^  at  both 
very  small  4  very  large  4- 

From  Equation  (19),  the  solution  for 

Fq(C)  =  erf(4)  (23) 

For  small  value  of  4>  approximated  as 


2  “  2 
^  '  7n  .=0  ^  ' 


Using  Equation  (24),  Gj(^)  was  determined  from  Equation  (21). 

G  ' (t)  =  -V2F  -  V24F 


■Jn  -Jn 


Thus , 

G-(4)  ~  -  2  (^)  ^  (26) 

^  ■  Vn 

By  substitution  of  Equations  (24)  and  (26)  into  Equation  (20),  we  have 

24Fj'  +  Fj"  -  zl  i  -  ^  (27) 

Vn 

By  setting  ^  =  0  in  Equation  (27),  and  using  the  boundary  condition 


Fj(4=0)  =  0,  we  obtain 

Fj"(4=0)  2  0 


Thus ,  we  have 


f,-(E=o)  -  - 
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Equation  (20)  is  true  only  for  very  small  t*  due  to  the  approximation 

used  in  Equation  (24).  For  very  small  t*  >/2/V^t*  >  >  2^2/n  -Jt* ,  thus 

Equation  (29)  could  be  approximated  by  neglecting  the  terms  of  0(Vt*) 
3/2 

and  0(t*  ).  Finally,  for  very  small  t*, 


V2 

3nln=o  - 


(30) 


For  very  large  t*,  the  ^  term  in  Equation  (2)  becomes  negligible, 
and  Equations  (1)  and  (2)  become  the  Blasius  problem.  Its  solution  was 
shown  by  Schlichting  (1968)  and  Watkins  (1975)  as 


l^ln=0  -  ^  0.i,695  as  t*  »  1  (31) 

Equation  (31)  is  assumed  to  be  the  expression  for  very  large  t*. 

Using  Equations  (30)  and  (31),  the  method  of  Rosezweig  was  applied 

0  f 

and  an  approximate  form  for  ^jq=0’  expressed  by  Equation  (32), 

was  obtained.  The  derivation  of  Equation  (32)  is  shown  in  Appendix  B. 


^o'  =  lnln=0  -  {-2f'^(*)t*}  +  f'(»)  erf  {V2t*f'(»)}  (32) 

where  f ' (“)  =  0.4695.  Calculated  results  of  f^'  at  several  points 
along  the  t*-axis  are  tabulated  in  Table  2.  Using  these  f^'  values. 
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Equation  (32)  is  plottpd  in  Figure  12.  As  implied  hy  assumption  (c)  of 
the  time  dependent  flow  model,  infinitely  large  surface  shear  stress  at 
t*  =  0  are  unavoidable.  As  time,  t*,  increases,  approaches  f'(*)  = 
0.4695  asymptotically. 

3.3  Stochastic  Averaging  Techniques 

The  simplified  time  dependent  flow  model  discussed  in  the  last 
section  cannot  be  directly  used  due  to  the  random  nature  ot  the  sweep 
motion.  In  order  to  account  for  the  spatial  and  temporal  dependence  of 
the  bursting  process,  a  stochastic  averaging  technique,  which  incorpo¬ 
rates  the  simplified  model,  was  developed  to  predict  the  probability 
density  distribution  of  the  surface  shear  stress  under  a  turbulent 
boundary  layer. 

This  stochastic  averaging  technique  starts  with  the  construction 
of  a  histogram  of  £q' •  A  general  criterion  for  the  truncation  of  Equa¬ 
tion  (32)  is  not  available,  so  truncational  values  of  t*  were  estimated 
empirically  by  utilizing  the  experimental  results  to  be  discussed  in 
Chapter  V.  By  expressing  t*  as  t'-'  =  t  U^./x,  truncational  values  of 
t*  could  be  estimated  by  considering  a  sweep  motion  which  originated 
at  the  location  of  the  instability.  According  to  the  proposed  flow 
model,  this  specific  sweep  motion  would  affect  the  flow  variables  near 
the  surface  until  the  next  sweep  motion  starts.  Thus  for  a  simplified 
case  for  which  the  frequency  of  occurrence  and  strength  of  the  sweep 
motion  are  constant,  the  period  between  two  sweep  motions  and  the 
strength,  U^,  together  with  x  should  be  used  to  obtain  the  truncational 
t*.  For  an  actual  flow  problem,  mean  bursting  period  T  and  most 
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As  concluded  from  the  study  of  Cliff  and  Sandborn  (1973),  the 
convective  velocity  measured  at  the  height  where  it  has  the  same  magni¬ 
tude  as  the  local  mean  velocity  was  considered  as  the  most  probable 
velocity  of  the  sweep  motion.  Results  summarized  by  Rao  et  al.  (1971) 
and  Falco  (1980)  for  TU^/6  observed  by  many  researchers  were  used  to 
estimate  T.  The  space-time  cross-correlation  measurements  of  velocity 
and  surface  shear,  which  are  discussed  in  Section  5.4  of  Chapter  V, 
provided  the  estimation  of  the  averaged  x-coordinate .  By  adopting 
U^T/6  =  7.5  and  using  the  results  shown  in  Table  6  for  6,  Figure  26 
for  (U  )  and  Figure  28  for  x,  the  following  values  were  obtained: 
a)  For  =  8.55  m/sec, 

6  =  6.528  cm  x  ~  6.0  cm,  (U  )  =4.3  m/sec, 

c  mp 

the  truncation  value  of  t*  is 

T(U  )  7.5  6  (U  ) 

t*  =  —piE  s - =  4.10 

X 


00 


b)  For  U  =  10.2  m/sec 


6  -  6.350  cm,  x  =  7.0  cm, 
the  truncation  value  of  t"  is 


(U  )  =5.1  m/sec , 

c  mp 


7.5  6  (U  ) 

t*  = - ^  .  3.41 

U  X 
00 


Based  on  the  above  estimations,  t*  =  4.0  was  chosen  as  a  representa¬ 
tive  truncation  time  for  applying  Equation  (32)  to  the  following 
demonstrat ion. 

Equation  (32)  between  f-'  =  0.0  and  t*  =  4.0  was  plotted  with 
200  equal  intervals  of  At*  =  0.02  on  the  t*-axis.^  For  each  fp'  a 

^Principles  of  digital  method  for  estimate  of  probability,  Bendat  and 
Piersol  (1971),  was  applied  in  the  present  study. 
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window  of  width  =0.1  was  used,  the  number  of  At''  intervals, 
including  fraction  of  an  interval,  which  has  a  value  of  f^'  within  a 
specific  window  was  estimated  graphically.  This  number  divided  by  the 
total  number  of  At*  interval  and  window  width,  200  and  0.1  respec¬ 
tively,  resulted  in  the  estimated  probability  density  for  f^'  covered 
in  a  specific  f^'  window.  The  above  procedure  is  expressed  mathemati¬ 
cally  as  follows; 
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where 
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=  width  of  each  f^'  window 
=  central  point  of  ith  f^'  window 
=  total  number  of  At*  interval 
=  200  in  the  present  study 
=  number  of  At*  interval  such  that 
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=  estimated  probability  density  for 
ith  fp'  window. 
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The  results  obtained  using  the  procedure  are  tabulated  in  Table  3,  and 

the  estimated  probability  density  function  is  plotted  in  Figure  13. 

The  relation  between  surface  shear  stress  T  and  f  „ '  ,  i.e., 

w  0 


T 

w 


MU, 
V2  X 


(34) 


was  obtained  by  introducing  x*^,  y°,  u*^  and  n  >  (which  were  defined  lu 
the  last  section)  into  the  relation  tor  f^' 


M.  _  ^u" 
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(35) 


31 


The  interrelation  among  T^,  and  x  was  unknown.  For  the  present 
study,  it  was  assumed  that  for  a  specific  boundary  layer  problem  the  x 
value  could  be  estimated  using  the  most  probable  T  and  U  ,  which 
were  determined  experimentally,  and  the  most  probable  f^'  value  cal¬ 
culated.  Consequently,  Equation  (34)  was  simplified  as 


(T  )  =  u)  (U  (f.-) 

w  mp  c  mp  0  mp 


(36) 


where 


u  1  ^  . 

lu  =  -  =  constant 

•J2  -Jxv 


(37) 


depends  on  a  specific  problem,  which  would  be  determined  from  the  most 
probable  values  of  and  for  a  specific  problem. 

By  assuming  that  the  probability  density  function  of  for  a 

specific  problem  is  known  and  that  axis  was  divided  into  J  windows 

of  equal  width  of  AU^,  the  probability  density  Pjt’',)^]  for 

withing  the  jth  window  could  be  estimated.  used  for  the 

grid  point  on  the  axis  where  its  probability  density 

to  be  calculated.  From  Equation  (36)  and  (37),  for  the  jth,  window, 

there  is  a  particular  ith  window  of  f^'  which  would  generate  a 

interval  covering  Thus  the  overall  probability  density 

due  to  the  above  statistical  U  distribution  could  be  estimated  by 
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(U  ) .  =  jth  U  window, 
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J  =  total  number  of  U  windows, 
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K  =  total  number  of  (T^)j^  points. 
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=  ith  window  of  f„'  which  generates  a  T  interval 
0  1  0  w 

covering  (t  ).  when  associated  with  (1)  ). 

w  k  t  j  , 

=  ml  (U^).  Af^-, 

(t  ).  =  grid  point  on  r  axis  at  which  the  probability 

calculated. 

A  crude  estimate  of  T  was  obtained  by 

w 


T  =  1  (t  ).  •  P  (u  )  •  A(t  ). 

w  ,  ,  w  k  w  k  w  k 

k=l 


The  flow  chart  for  the  above  procedure  of  probability  density 


calculation  for  is  presented  in  Appendix  C. 

To  demonstrate  the  above  stochastic  averaging  technique  and  its 
application  together  with  the  flow  model  to  predict  the  surface  shear 
stress,  several  sets  of  computations  were  made.  Experimental  results 
for  the  surface  shear  stress  for  the  cases  of  ll  =  10.2  m/sec  and  U„ 

00  OO 

=  8.55  m/sec,  (which  are  shown  in  Figures  23  and  26),  were  used  to 


obtain  the  required  constant  w.  The  convective  velocity  measured  at 
the  height  where  the  convective  velocity  and  the  mean  velocity  were 
equal  was  used  as  the  most  probable  velocity  of  the  sweep  motion.  The 
following  values  were  estimated  from  the  experimental  data: 

a)  For  U  =  8.55  m/sec, 

(U  )  =  A. 30  m/sec,  (x  )  =  0.1  N/m^ 

c  mp  ’  w  mp 

b)  For  U  =  10.20  m/sec, 

oo 

(U  )  =  5.10  m/sec,  (t  )  =  0.153  N/m^. 

c  mp  w^mp 

By  assuming  “  0.5,  U)  =  0.02666  and  0.0224  were  obtained  for  the 

above  two  cases  respectively. 

Using  these  values  and  assuming  several  different  possible 
distributions  for  the  sweep  velocity,  results  as  shown  in  Figures 


iJ 


14a)-l4i)  were  obtained.  The  corresponding  experimental  measurements 
are  also  shown  in  these  figures  for  comparison.  A  Gaussian  probability 
distribution 


_  ,  (U  -  (U  )  I*" 

P(U  )  =  exp  [ - ^ - - -  ) 

20 


and  a  modified  Rayleigh  distribution, 


,  “c  ■  ,  I'Jc  -  ‘"c^rnlnl  , 

P(U,)  = - 2 - exp  1 - -1 - 1 

c  2c 


were  employed  in  the  present  study.  In  the  above  equations,  ^'^c^mp 
was  used  as  the  mean  value  for  the  Gaussian  distributions,  (13  )  .  and 
c  were  selected  such  that  the  assumed  Rayleigh  distributions  would 
have  the  most  probable  approximately  equal  to  the  '^hich 
corresponds  to  determined  from  Equation  (36). 

Computed  results  show  that  the  skewed  shape  of  the  probability 


density  distribution  of  surface  shear  stress  under  a  turbulent  boundary 
layer  is  predicted  for  each  case  tested.  The  predicted  magnitude  of 
the  most  probable  shear  stress  is  in  reasonable  agreement  with  the 
measurements.  Compared  with  the  measured  probability  distribution  of 
surface  shear,  the  analysis  over  predicts  the  high  magnitude  fluctua¬ 


tions.  This  over-prediction  is  due  in  part  to  the  singular  point 
encountered  at  t*  =  0  when  Equation  (32)  was  applied.  The  singular¬ 


ity  combined  with  the  summation  scheme  for  computing  Equation 
(39),  will  contribute  to  the  over  prediction  of  T^.  It  is  also  pos¬ 
sible  that  the  frequency  response  of  the  measuring  instruments  may  make 
the  results  at  the  large  magnitudes  questionable. 


It  was  possible  to  improve  the  agreement  (less  than  ±  20%)  between 
the  measured  and  predicted  values  of  surface  shear  by  employing 


arbitrary  values  of  the  standard  deviation  and  the  equivalent  velocity 
distribution. 

Increase  of  the  standard  deviation  for  the  Gaussian  distribution 
and  its  equivalent  quantity,  c,  for  the  modified  Rayleigh  distribution 
reduces  the  magnitude  of  the  peak  and  flats  out  the  shape  of  the  re¬ 
sulting  probability  curve.  The  flatting  effect  is  more  apparent  toward 
the  lower  shear  stress  side  than  toward  the  higher  side,  as  can  be  seen 
by  comparing  Figures  I4h)  and  c)  with  Figure  lAa),  and  Figures  I4e)  and 
f)  with  14d).  These  results  seems  to  justify  the  approach,  which  leads 
to  Kquation  (32J,  as  well  as  the  proposed  physical  modc'l  which  predicts 
a  Rayleigh  type  surface  shear  stress  response  to  the  bursting  motion. 
Moreover,  the  agreement  between  the  measured  and  predicted  values, 
shown  in  Figures  I4a)  and  l4d) ,  suggests  that  the  distribution  of  sweep 
velocity  could  be  similar  to  a  Gaussian  distribution  with  ratio  of  mean 
value  to  standard  deviation  of  .about  0.10. 

Parameters  used  in  these  computations  and  the  mean  surface  shear 
stress  obtained  are  summarized  in  Table  4. 


CHAPTER  IV 


EXPERIMENTAL  SETUP  AND  PROCEDURE 

The  experimental  study  was  performed  at  the  Engineering  Research 
Center  of  Colorado  State  University.  Measurements  were  made  in  a  small, 
variable  geometry,  wind  tunnel  in  a  region  where  the  longitudinal 
pressure  gradient  was  nearly  zero. 

The  major  part  of  the  measurements  consist  of  correlations  of 
surface  shear  stress  and  fluctuating  velocity  in  the  mean  flow  direc¬ 
tion.  Convective  velocity  and  probability  density  of  turbulent  velo¬ 
city  and  surface  shear  stress  were  also  evaluated.  The  signals  from  a 
dual  wire  probe  and  a  surface  hot  wire  were  recorded  using  a  FM  re¬ 
corder.  A  correlation  and  probability  analyzer  was  used  to  evaluate 
the  cross-correlation  between  each  set  of  signals  and  the  probability 
distribution  for  the  individual  signals.  The  cross-correlation  and 
probability  of  the  voltage  signals  were  plotted  directly  using  an  X-Y 
plotter.  The  instrumentation  employed  for  these  measurements  included 
constant-temperature ,  hot-wire  anemometers,  an  analog  spectrum  ana¬ 
lyzer,  filters  ,  an  analog  percentage  time  analyzer,  and  a  HP-1000 
computer  facility.  Pressure  and  velocity  measurements  were  evaluated 
using  a  capacitance  pressure  transducer. 

4.1  Wind  Tunnel 

An  open-return,  vitriable  geometry,  wind  tunnel,  shown  in  Figure  15 
was  used.  The  test  section  was  made  of  plexiglass  and  has  a  cross 
section  of  45  cm  x  45  cm.  A  1.37-m  fan  driven  by  a  variable  speed 
motor  was  located  downstream  of  the  test  section.  The  free  stream 
velocity  was  controlled  manually.  Time  variations  in  the  free  stream 
pressure  were  limited  to  ±3  percent. 
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The  test  section  consists  of  a  zero  pressure  gradient  region  of 

approximately  190  cm  in  length  over  which  a  turbulent  boundary  layer  of 

the  order  of  6  cm  in  thickness  was  developed.  The  equivalent  momentum 

thickness  Reynolds  number,  R.,  is  of  the  order  of  4.0  x  10^.  This  flat 

0 

plate  region  was  followed  by  a  curved,  adverse  pressure  gradient 
region.  Figure  15  shows  the  setup  of  the  test  sectioii. 

The  flow  measurements  were  limited  to  the  zero  pressure  gradient 
region.  At  the  centerline  of  the  wind  tunnel  and  198  cm  downstream  of 
the  entrance  of  the  test  section,  a  surface  hot  wire  was  m(.>uuLed  flush 
to  the  formica  surface.  A  series  oi  static  pressure  holes  of  0.05  cm 
diameter  were  drilled  along  the  center-line  in  order  to  evaluate  the 
static  pres.sure  distribution.  An  actuator  monnicnl  beneath  tlie  ()!ati 
was  used  to  traverse  the  probes  through  the  bountlary  l.ayc't 

The  coordinate  system  used  is  that  the*  x-axi:.  is  patille!  to  the 
surface  and  the  flow  (positive  x  is  measurc’d  upstream  ot  ‘.he  surface’ 
hot-wire),  the  y-axis  is  perpendicular  to  the  surfne  and  the  flow,  ami 
the  z-axis  is  parallel  to  the  surface  and  per|iei  >1  i  i  u  I  a  r  to  the  flow 
direction . 

4.2  Instrumentation 

Hot  wi  re_protK’.s .  Hot  wire.'--  weie  used  as  the  sensing  element  for 
all  correlation  measurements,  probability  me.i  sii  rement  s  and  surface 
shear  stress  evaluation.  The  material  lor  the  hot  wires  was 
pi  at  i  num-8‘’4  tungsten,  t).OI  mm  diameter  for  all  sensors.  The  length  of 
wires  used  ranged  t  roiii  0,03  cm  to  O.Oh  cm.  The  wires  were'  operated  by 
constant  temperature,  hot  wire  anemomc'tf'r  ciiiuils. 

Hot  wire  sensors  of  three  different  designs  were  employed.  They 
were  single  wire  probes,  a  dual  wire  probe  and  a  surface’  wire.  The 


single  wire  probe  was  used  tor  the  turbulent  measurements  since  it 
could  be  traversed  very  close  to  the  wall.  The  dual  wire  probe,  which 
had  the  approximate  dimensions  shown  in  Figure  16a),  was  used  mainly 
for  the  measurements  of  convective  velocity,  and  the  cross  correlation 
of  surface  shear  stress  and  fluctuating  velocities  through  the  boundary 
layer.  The  rear  wire  signal  of  the  dual  probe  was  correlated  with  the 
lead  wire  signal  in  order  to  measure  a  convection  time  for  the  turbu¬ 
lent  structure  to  travel  from  the  lead  to  the  rear  wire.  The 
information  between  the  surface  shear  stress  and  related  motion  in  the 
boundary  layer  was  obtained  by  measuring  the  cross-correlation  between 
the  probe  hot  wires  and  the  surface  hot  wire  signals. 

The  surface  hot  wire  was  used  to  measure  the  fluctuating  surface 
shear  stress.  It  was  constructed  by  mounting  the  wire  directly  on  the 
surface,  such  that  the  heat  transfer  from  the  wire  was  limited  to  the 
linear  velocity  region  of  the  boundary  layer. 

Correlation  and  probability  analyzer.  For  correlation  and 
probability  measurements,  a  Signal  Analysis  Industries  Corporation 
Correlation  and  Probability  Analyzer,  Model  SAI-42  was  used.  The 
correlation  analyzer  was  a  hybrid  computer,  which  uses  both  analog  and 
digital  techniques.  Correlation  analysis  provides  a  quantitative 
measure  of  the  degree  of  similarity  between  signals  as  they  appear 
relative  to  one  another  in  time.  The  SAI-A2  model  provides  auto-  and 
cross-correlation  functions  with  incremental  lag  or  time  delays  values 
ranging  from  1  psec  to  1  sec  resulting  in  total  time  delays  of  from  100 
psec  to  100  sec.  An  auto-  or  cross-correlation  function  was  determined 
simultaneous ly  at  100  incremental  lag  points  so  that  a  complete  corre¬ 
lation  function  is  displayed  at  one  time. 
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The  probability  analysis  provided  information  concerning  the 
likelihood  that  the  amplitude  of  a  signal  lies  within  a  prescribed 
interval.  The  analyzer  determined  probability  at  100  continuous 
amplitude  intervals.  Value  ol  each  amplitude  bound  w.is  obtained 
through  a  calibration  procedure,  which  employed  a  sine  wave  of  known 
ampl itude . 

Tape  tran^spo^.  A  fre<piency  modti  1 .1 1  i  o(i  type  o(  tape  recorder  was 
used  for  the  study.  Seven  channels  could  he  used  to  record  .signals 

simultaneously,  however,  only  the  o(id  iiuiiihet  were  used  t  >  avoid  the 
possible  time-displacement  errors  cmised  hy  the  idiysiial  separation  of 
the  odd  and  even  numbered  heads. 

Time  percentage  analyzer.  A  analog  lime  p{'rcetil.age  analyzer  which 
was  designed  hy  Finn  and  Sandborn  (196A)  was  used  to  measure  the  inlor- 
mittency  through  the  boundary  layer.  This  analyzer  measures  tlu' 
percent  of  time  an  ac  signal  is  non-zero. 

^  pitot  tube  with  an  eliiptical  shaped 

opening  and  the  dimensions  shown  on  Figure  16h)  was  used  to  measure  the 
local  mean  velocity.  The  output  of  the  pilot  tube  pressure  was  applied 
to  the  high  side  of  a  diilerential  pressure  transducer,  while  tiie  low 
side  was  connected  to  the  local  surlacre  static  pressure  tap.  The 

pressure  was  sampled  and  digitized  at  a  rate  of  600  Hz.  and  averaged  for 
20  sec  by  d  HP-IOUO  computer,  Each  dynamic  pressure  reading  was  con¬ 
verted  to  the  equivalent  velocity  and  then  averaged.  The  averaged 
velocity  and  root -mean-squa re  value  of  the  lliictuating  component  were 
returned  to  the  keyboard  terminal  instantaneously.  Using  this  small 
size  pitot  tube,  the  mean  velocity  at  a  hc'ight  of  y""  >  2.0  above  the 


wall  w.is  determined. 


Height  indicating  system.  To  explore  the  flow  properties  in  the 
sublayer  region,  the  aitual  height  of  the  measuring  probe  must  be 
accurately  determined.  File  measuring  probe  was  mounted  on  a  rod,  which 
was  connected  to  a  motor-driven  aitu.alor.  Attached  to  the  actuating 
mechanism  was  a  dial  indicatoi  v'lh  .m  aciur.ity  of  0.02.^  cm  (0.001 
in.).  This  indicator  was  used  to  indie. i  -  the  height  of  the  measuring 
probe  above  the  surface.  A  thin  graphite  film  (pencil  mark)  was  ap¬ 
plied  on  the  flat  plate.  .An  electrical  circuit  was  used  to  indicate 
when  the  probe  was  in  contact  with  the  graphite  film,  so  the  exact 
location  of  it  could  be  determined.  The  dial  indicator  was  then  set  to 
tne  effective  height  of  the  sensor  (hot  wire  or  pitot  tube)  above  the 
surface.  A  special  jig  with  a  calibrated  optical  microscope  was 
employed  to  determine  the  probe  dimensions.  Correction  for  velocity 
gradient  and  averaging  errors  were  made. 

4.3  Evaluation  of  the  Hot-wire  Signals 

Time  series  of  signals  from  the  dual  hot-wire  and  the  surface 
hot-wire  were  recorded  simultaneously  for  3  minutes.  A  sine  wave, 
which  was  used  both  as  identification  and  for  calibration  was  recorded 
before  each  set  of  measurements.  The  hot-wire  voltages  were  controlled 
within  specific  limits  to  insure  that  the  recordings  contained  all  the 
frequency  information  from  DC  to  the  upper  limit  of  the  sensors.  The 
sine  wave  was  employed  to  calibrate  the  gains  of  the  record  and  repro- 
iluce  amplifiers  of  the  tape  recorder. 

For  convective  velocity  evaluation,  the  longitudinal  separation 
between  the  lead  and  the  back  wires  were  divided  by  the  time  del.iy  of 
the  maximum  correlation.  For  a  detailed  discussion  of  the  evaluation 
of  convective  velocities,  the  reader  is  referred  to  Cliff  and  Sandborn 


CHAPTER  V 


EXPERIMENTAL  RESULT  AND  DISCUSSION 
Measurements  were  made  at  free  stream  velocities  of  approximately 
10.20  m/sec,  9.36  m/sec  and  8.55  m/sec.  The  corresponding  Rg  values 
are  ^050,  3500  and  3080,  respectively.  The  primary  part  of  the  mea¬ 
surements  was  related  to  the  surface  shear  stress  evaluation  at  the 
center  line,  196.6  cm  downstream  of  the  entrance  of  the  test  section. 
The  location  (196.6  cm)  was  chosen  as  the  origin  of  the  coordinate 
system.  Also,  the  positive  streamwise  coordinate  (x-coordinate)  was 
taken  as  upstream  of  the  196.6  cm  station  in  the  wind  tunnel.  At  Rg  = 
3080,  modification  of  surface  shear  stress  by  blocking  the  passage  of 
the  sweep  motion  was  investigated.  This  investigation  was  aimed  to 
justify  the  proposed  flow  model. 

5.1  The  Flow  Field  over  the  Test  Model 

The  time  averaged  variables  of  the  turbulent  boundary  layer 
developed  over  the  test  section  were  evaluated.  The  object  was  to 
assure  a  nearly  zero  pressure  gradient  condition,  a  uniform  flow  and  a 
normally  behaved  boundary  layer  over  the  surface. 

The  static  pressure  distribution  along  the  centerline  of  the  test 
model  was  obtained  by  measuring  the  pressure  output  from  the  0.05  mm 
diameter  holes  referenced  to  the  value  at  a  location  of  120  cm  from  the 
entrance  of  the  test  section.  The  results  are  tabulated  in  Table  5  and 
are  plotted  in  Figure  17.  The  pressure  is  nearly  constant  from  30.0  cm 
to  200  im.  Although  there  is  a  slight  variation  in  pressure  along  the 
surface,  the  deviations  are  of  the  same  order  as  the  accuracy  that 
individual  reading  can  be  made. 
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At  a  tree-stream  veloeily  ot  -  12  80  m/sei  ,  mean  velocitv 

profiles  were  roeasureil  for  several  lateral  locations  at  x  =  -10.2  cm 
-15.20  cm  and  -24.11  cm,  as  shown  in  Fignre  18.  Velocity  variations  in 
the  lateral  direction  across  the  test  section  were  within  approximately 
±  1%.  The  profiles  measured  at  x  =  -24.11  cm  indicate  an  acceleration 
which  corresponds  to  a  favorable  pressure  grailieiiL  as  indicate!  on 
Figure  17.  The  flow  field  up  to  the  measuring  point  has  a  near  x.ero 
pressure  gradient  and  it  is  uniform  across  the  span. 

At  the  location  of  x  =  0  ami  z  =  0,  detailed  time  ave''.igei1 
properties  including  mean  velocity  profiles,  longitudinal  turbulf'u'o 
intensity  profiles  and  the  internn  1 1  ency  were  measured.  Talile  6  lists 
the  related  boundary  layer  (juantilies  calculated  from  the  mean  velocity 
measurements.  Also  listed  are  the  skin  friclion  rimf  f  iejent ,  ,  t/ie 

corresponding  surface  shear  stress,  i  ,  and  shear  stress  velocity,  II  . 
The  skin  friction  coefficients  were  calculated  from  the  one  parameter 
equation  of  Bell  (197d),  shown  in  Figure  19.  The  inner  region  similar¬ 
ity  plot  for  mean  velocity  is  shown  i  t\  Figure  20,  and  is  compared  with 
a  typical  "Law  of  the  w,i  I  1  .  ”  The  mean  velocity  profiles,  turbulence 
intensity  profiles  and  inlermi ttency  profiles  are  summarized  in  Figure 
21.  The  related  data  for  mean  velocity,  turbulence  intensity  and 
int ermittency  are  tabul.ited  in  Tables  7,  8  and  9,  respectively.  In 
Figure  2lh)  turhulenee  measurement  .s  w.ide  by  KIc’banoff  and  by  Kckelmann 
in  the  wall  region  over  flat  plates  arc  also  compared  with  the  present 
data . 


As  shown  in  Figure  20,  the  comparison  of  results  of  the  present 
study  with  the  law  of  the  wall  indicates  that  the  flow  corresponds 
closely  to  a  flat  plate  flow.  The  turbulence  intensity  measurements  of 
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the  present  study  is  bracketed  by  the  measurements  of  Klebanoff  and 
Eckelmann,  as  illustrated  on  the  insert  of  Figure  21b).  The  decrease 
of  turbulence  intensity  toward  the  wall  in  sublayer  region  0  <  yU^/v  < 
5  indicates  that  high-speed  fluid  does  reach  the  wall  and  travel 
adjacent  to  the  wall  as  suggested  by  the  time  dependent  model  for  the 
surface  shear  stress.  Figure  21c)  shows  a  comparison  of  intermittency 
measurements  of  the  present  study  with  the  relation  suggested  by 
Klebanoff.  Figure  22  shows  spectrum  measurements  at  several  heights 
across  the  boundary  layer.  These  plots  show  that  the  contribution  of 
eddies  in  the  higher  frequency  range  increase.s  as  the  wall  is 
approached . 

5.2  The  Probability  Distributions  of  Hot-wire  and  Surface-wire  Signals 
The  probability  density  distribution  of  the  surface  shear  stress, 
P(t^,),  and  the  streamwise  turbulent  velocity  components,  P(u'),  at 
several  heights  from  y*  ~  2.0  up  to  y*  ~  800  were  evaluated.  The 
probability  density  distribution  of  a  specific  quantity  provides  infor¬ 
mation  on  the  magnitude  of  the  fluctuating  values.  It  would  be  ex¬ 
pected  that  the  velocity  probability,  P(u'),  measured  in  the  sublayer 
region  should  have  a  similar  pattern  as  the  surface  shear  stress, 
P(T^)  .  l.arge  magnitude  excursions  of  the  shear  stress  and  the  stream- 
wise  velocity  flucutations  have  been  reported  near  tlie  wall. 

Due  to  the  non-linear  relationship  between  the  magnitude  of 
(  surface  shear  stress  and  transducer  output  (voltage),  a  correction  is 

j  required.  The  technique  for  evaluating  the  surface  shear  stress  with 

I  surface  hot  wires  is  discussed  by  Sandborn  (1979).  Figure  23  shows  the 

I 

'  corrected  curves  of  P(T  )  for  =  8.55  m/sec  and  U  =  10.2  m/sec. 

'  ^  y '  00  00 

The  surface  shear  stress  is  skewed  with  large  positive  magnitude  fluc- 
j,  tuation,  as  great  as  2.75 
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The  probability  density  distribution  of  the  streaniwise  velocity 
fluctuations  at  various  heights  are  shown  in  Figure  24.  The  abscissa 
of  these  plots  was  normalized  with  the  local  mean  velocity  U.  The 
variations  of  the  velocity  probability  distribution  with  height  are  in 
agreement  with  similar  results  of  Brodkay,  Wallace  and  Kckelmann 
(1973).  In  the  sublayer  region,  (y""  <  10),  P(u')  is  skewed  with  large 
positive  values.  The  maximum  fluctuation  observed  at  y*  ~  2.0  is  n( 
the  same  order  of  magnitude  as  the  local  mean  velocity  F.  The  prob¬ 
ability  distributions  approach  a  Gaussian  sliape  and  gradually  ski-w  to 
the  opposite  direction  with  vertical  height,  y'  >  130,  and  thf 
distribution  has  a  very  small  standani  deviation. 

The  negative  side  and  positive  side  of  P(u')  could  he  viewed  as 
representing  flow  of  speed  lower  or  higher  than  the  mean  velocity 
respectively.  For  y'"  <  10,  flow  of  speed  higher  than  the  mean  velocity 
occurs  over  a  small  percentage  of  the  total  time.  The  high-speed  flow 
is  then  gradually  retarded.  The  retarding  process  lasts  longer  than 
the  acceleration.  Brodkey,  Wallace  and  Eckelmarin  (1974)  were  able  to 
determine  the  individual  contribution  of  each  motion  in  a  bursting 
cycle  to  the  retarding  and  accelerating  process.  Their  results  for  y'-' 
=  3.4  is  illustrated  in  Figure  25.  it  shows  that  both  the  sweep  motion 
and  the  outward  interaction  are  accounted  for  in  the  tlow  acceleration; 
hut  the  sweep  is  more  effective  than  tlie  outward  interaction  and  is 
responsible  for  the  much  higher  magnitude  velocity.  It  also  indicates 
that  the  ejection  motion  contributes  more  to  the  deceleration  of  (hi' 
flow  than  the  inward  interaction  does. 

Comparison  of  the  shape  of  P(t^)  and  P(n' )  in  the  sublayer  region 
reveals  a  similar  pattern.  This  similarity  demonstrates  the  interrela¬ 
tion  between  the  sublayer  motion  and  surface  shear  stress.  Although 


the  maximum  fluctuation  of  u  is  only  of  the  order  of  2U,  maximum 
fluctuations  of  the  surface  shear  stress  approache  3  T  .  This  increase 
in  magnitude  demonstrates  the  importance  of  the  non-linear  response 
involved  in  a  turbulent  boundary  layer. 

5.3  Convective  Velocity  Measurements 

The  dual-wire  probe  was  used  for  the  measurements  of  convective 
velocity  over  the  test  model.  The  two  wires  were  separated  by  a  dis¬ 
tance,  AW  which  ranged  from  0.550  cm  to  0.625  cm.  Convective  velocity 
profiles  measured  at  x  =  0.0  cm,  z  -  0.0  cm  for  =  10.20  m/sec  and 
8.55  m/sec  are  shown  in  Figure  26  together  with  the  corresponding  mean 
velocity  profiles. 

The  convective  velocity  was  smaller  than  the  mean-flow  velocity  in 
the  outer  region.  But,  the  convective  velocity  was  much  higher  than 
the  mean  velocity  in  the  sublayer  region.  Since  convective  velocity  is 
a  measure  of  the  rate  at  which  the  turbulent  structures  are  convected, 
results  shown  in  Figure  26  imply  that  turbulent  structures  in  the 
sublayer  are  transported  at  a  rate  faster  than  the  local  mean  velocity. 
This  is  possible  only  when  a  high-speed  motion  exists  in  the  sublayer 
region.  Measurements  of  convective  velocity  thus  provides  direct 
evidence  for  the  existence  of  high-speed  velocity  associated  with  the 
sweep  motion  in  the  sublayer  region,  as  proposed  in  the  physical  model. 

5- A  Cross-correlation  between  Surface  Shear  Stress  and  Streamwise 

Turbulent  Velocity 

Space  and  time  correlations  of  the  time  series  of  boundary  layer- 
hot-wire  and  surface-hot-wire  signals  were  evaluated.  At  selected 
heights,  the  surface-hot-wire  and  boundary-layer-hot-wire  signals  from 
the  du.ii  wires  were  recorded  simultaneously  for  at  least  3  minutes. 
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The  cross-correlal ion  between  the  surface  shear  stress,  T  ,  and  the 

w 

streamwise  turbulence  velocity,  u,  defined  as 


R  (Ax,  Ay,  Ax;  riAt )  =  1  im  T  (t)  u(t  +  nAt)dt  (42) 

V  T  -0-  s  0 

s 

was  evaluated  using  the  hybrid  analog  digital  correlator.  In  the  above 
equation.  Ax,  Ay  and  Az  represent  separation  of  the  hot  wire  and  sur¬ 
face  wire  in  the  x,  y  and  z  directions,  respectively.  nAt  is  the  time 
increment  and  is  the  total  summation  time  used  in  tiie  computation  oi 

the  correlation.  R  for  n  =  I  to  n  =  100  were  calculated.  in  liiis 

T  u 
w 

study.  At  ranging  from  20  x  10  sec  to  500  x  10  ^  sec  and  ranging 
from  30  sec  to  66  sec  were  used. 

The  major  portion  of  the  correlation  measurements  were  inadi’  iiy 
aligning  the  hot  wire  with  the  surface  wire  in  the  streamwise  direction 
and  keeping  z  =  0.  The  location  of  the  hot  wire  varie.s  from  x  =  9.0  cm 
to  X  =  -7.0  cm.  At  each  location,  measurements  were  made  at  several 
heights  from  y  ~  0.008  cm  up  to  y  =  1.27  cm.  The  equivalent  non- 
dimensional  heights  ranged  from  y''-  =  3  to  Y"  -  250.  For  the  z  =  0 
rase,  a  single  maximum  of  the  curve,  (K  )  is  nearly  a  constant 


w 


value  for  a  range  of  nAt,  thus,  it  was  ilifficult  to  determine  a  unique 


Value  of  At  This  uruertaintv  in  At  increases  as  llie  magnitude 

max  •  max 

of  (R  )  becomes  smaller.  In  the  present  study,  the  midpoint  of 

T  u  max 
w 

the  uncertain  range  of  nAt  was  chosen  as  (At)  .  The  data  of  (R,  ) 

”  max  T  u  max 

w 


and  (At)  are  presented  in  Appendix  D. 


(R  )  and  (At)  obtained  at  each  measuring  point  in  the 

T  u  max  max 

w 


domain  about  the  surface  shear  stress  sensor  were  employed  to  construct 
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the  peak  isororrelat ion  curves  and  the  peak  time  delay  curves  shown  on 


Figures  28  and  29,  respectively. 


Three  possible  errors  were  investigated  in  the  process  of  the 


correlation  measurements.  One  is  due  to  the  unsteadiness  of  the  flow, 


another  is  related  to  the  repeatability  of  the  results  and  the  third  is 


the  phase  miss-match  between  the  measuring  instruments.  Evaluation  of 


(R,  )  for  different  segments  of  the  time  series  of  the  hot  wire 

T  u  max 
w 

and  the  surface  wire  signals  for  a  specific  value  of  Ax  and  Ay  separa¬ 


tions  showed  that  the  maximum  deviation  of  individual  measurements  from 


the  mean  value  of  (R,  )  is  ±6%  of  the  mean.  To  check  the  repeat- 

T  u  max 
w 

ability  of  the  results,  comparison  of  results  obtained  at  the  same  lo¬ 


cations  on  different  days  were  made.  On  one  day,  the  measurements  were 


made  using  the  lead  wire  of  the  dual  wire  probe;  and  on  the  second  day. 


the  back  wire  of  the  dual  wire  probe  was  used.  The  comparison  is  il¬ 


lustrated  in  Figure  30.  From  this  figure,  a  maximum  deviation  of 


day-to-day  measurements  of  less  than  ±10%  was  estimated.  Based  on 


these  observations,  the  possible  maximum  uncertainty  related  to  an 


individual  evaluation  is  less  than  ±16%  if  the  two  uncertainties  add 


linearly. 


As  illustrated  in  Figure  28,  the  peak  isocorrelat i on  curves  for 


the  three  R^  values  show  common  characteristics.  The  isocorrelation 

lines  for  (R  )  =  0.40,  0.50,  0.60,  0.70  and  0.80  form  closed  loops 

w 

around  the  surface  sensor.  A  region  of  relatively  low  (R  )  ,  0.2  < 

"  X  u  max 

w 

(R  )  <0.3,  which  is  designated  by  the  symbol  B  exists  downstream 

w 

of  the  contour  of  (R  )  =0.40. 

X  u  max 
w 

Before  discussing  the  physical  meaning  of  the  present  peak 


isocorrelation  curves,  the  correlations  between  surface  shear  stress 
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and  turbulent  velocity  in  a  duct  evaluation  of  Rajagopalan  and  Antonia 
(1980)  are  of  specific  interest.  A  set  of  results  obtained  by 
Rajagopalan  and  Antonia  are  shown  in  Figure  31,  in  which  the  subscripts 
1  and  h  represent  the  low-pass  filtered  frequencies  and  the  high- 
pass  filtered  frequencies  of  a  signal.  These  plots  illustral<'  that  tor 
the  turbulent  velocity  u  measured  outside  the  inner  region  only  u^ 
correlated  with  T  ,  .  For  u  measured  in  the  inner  region,  the  cor- 
relation  of  and  u^  was  much  higher  than  the  other  frequency 
correlations . 

Based  on  the  above  finding  of  Rajagopalan  and  Antonia,  we  can 
interpret  the  isocorrelation  curves  shown  in  Figure  28  as  follows;  in 
the  outer  region  the  correlations  represent  the  large  scale  motion 
related  to  the  surface  shear  stress;  while  in  the  inner  region,  the 
correlation  between  the  entire  u  and  signals  are  present,  but 
the  correlation  between  the  high  frequency  part  of  u  and  is 
predominant.  The  high  frequency  part  of  the  turbulent  velocity  in  the 
inner  region  was  confirmed  to  be  associated  with  the  bursting  phenome¬ 
non  by  Kim  et  al.  (1971)  and  others.  Brown  and  Thomas  (1977)  observed 
that  the  large  magnitude  surface  shear  stress  is  also  at  the  high 
frequencies.  The  close  relation  between  the  sweep  motion  and  surface 
shear  stress  is  evident. 

The  peak  isocorrelation  curves  of  Figure  28  for  the  inner  region, 
y’’*’  <  30,  appear  to  be  nearly  parallel  to  each  other  between  x  ~  -  1.5 
and  X  ~  3.0.  On  the  downstream  side  of  the  surface  wire,  the  isocorre¬ 
lation  lines  for  .5  <(R  )  <  .6  are  distorted  to  an  elongated 

X  u  max 
w 

loo])  with  its  center  inclined  to  the  wall  at  a  small  angle.  ,Sincp  the 
high  frequency  part  of  u,  which  is  directly  related  to  the  bursting 


phenomenon,  is  the  major  contribution  to  the  correlate  with  the  surface 
shear  stress,  the  isocorrelation  lines  in  the  inner  region  must  indi¬ 
cate  the  approximate  passage  and  evolution  of  the  sweep  motion. 

The  peak  time  delay  curves,  shown  in  Figure  29,  provide 
information  about  the  evolution  of  the  large  scale  motion  in  the  outer 
region.  These  curves  start  with  a  nearly  vertical  but  irregular  shape 
in  the  upstream  region  and  gradually  incline  toward  the  wall  for  the 

downstream  direction.  The  peak  time  delay  curve,  (At)  =  0.0  sec, 

max 

represents  the  mean  locus  of  the  upstream  and  downstream  sides  of  the 
large  scale  motion,  as  discussed  by  Falco  (1978).  The  angle  between 
the  mean  of  (At)  =  0  curve  and  the  surface  is  13.0°,  11.13°  and 
12.53°  for  Rq  =  4050,  3500  and  3080  respectively.  These  values  are 
lower  than  the  value  of  approximately  18°  reported  by  Brown  and  Thomas 
(1977)  for  a  turbulent  boundary  layer  at  a  higher  Reynolds  number,  Rg  = 

1.06  X  10^. 

A  limited  number  of  measurements  of  R  (Ax,  Ay,  Az;  nAt)  with  Az 

w 

f  0  were  also  evaluated  for  Rg  =  4050.  The  major  part  of  these 
measurements  were  made  in  a  domain  of  x  from  0.0  cm  to  5.0  cm  and  z 
from  0.0  cm  to  1.5  cm  on  one  side  of  the  center  line  only.  The  hot¬ 
wire  probe  was  placed  at  y  =  0.127  cm,  0.254  cm  and  0.381  cm  above  the 
surface.  As  the  probe  was  moved  away  from  the  center-line,  the  magni¬ 
tude  of  correlation  decreased  rapidly,  however,  a  multimaximum  correla¬ 
tion  curve  was  obtained.  The  magnitude  of  correlation  became  negli¬ 
gible  for  z-distances  greater  than  1.5  cm. 

Figure  32  shows  typical  cross  correlation  curves  obtained  for  Ax  = 
5.0  cm.  Ay  =  0.127  cm  and  several  z  locations.  The  value  of  the 
maxima  measured  at  each  location  is  plotted  on  Figure  33  for  y  =  0.127 
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cm,  0.254  cm  and  0.381  cm.  For  Figure  33  the  data  with  solid  circles 
indicates  the  existence  of  a  multi-maximum  correlation  curve,  while 
those  with  open  circles  represents  the  existence  of  a  single  maximum 
correlation  curve.  An  approximate  boundary  between  these  two  groups  of 
data  points  is  constructed  for  each  plot  of  Figure  33.  This  boundary 
is  related  to  the  passage  of  the  large  scale  mol  ion.  The  lateral 
expanse  of  this  boundary  remains  nearly  constant  of  the  order  of  z'-  = 
100  from  y*  =  22.82  to  y*  =  84.7.  The  boundary  indicates  the  approxi¬ 
mate  scale  of  the  large  scale  motion. 

5.5  Surface  Shear  Stress  Modification 

Based  on  the  physical  model  discussed  in  Chapter  ill,  the  large 
magnitude  fluctuations  of  surface  shear  stress  are  directly  related  to 
the  sweep  motion  of  the  bursting  phenomenon.  Correlation  measurements, 
as  discussed  previously  in  section  5.4  of  the  present  chapter,  demon¬ 
strated  the  passage  of  the  sweep  motion.  It  was  postulated  that  if  the 
sweep  motion  could  be  modified  or  prevented  from  penetrating  directly 
into  the  sublayer,  the  large  magnitude  fluctuations  of  surface  shear 
stress  would  be  damped  or  eliminated. 

As  shown  it!  Figure  28c),  at  =  3080  the  sweep  effect  is  first 
evident  in  a  region  at  a  nondimensional  x-distrance  of  x*  ~  750  to 
1000,  and  at  nondimensional  y-distance  of  y*  ~  75  to  100.  The  random, 
three-dimensional  production  of  the  fluctuations  over  the  surface  may 
limit  somewhat  the  selection  of  a  specific  location  of  the  start  of  the 
sweep  motion.  If  the  instability  of  this  region  ran  be  blocked  from 
reaching  the  surface,  large  amplitude  surface  shear  fluctuations  should 
be  reduced  dramatically.  Thin  aluminum  plates  of  0.013  cm  thick  but 
with  different  widths  were  placed  at  different  heights  above  the 


surface  and  at  different  locations  with  respect  to  the  surface-hot-wire 
to  examine  their  effects.  Figure  3A  show  a  sketch  of  a  typical  plate. 
The  plates  have  a  marked  effect  on  the  magnitude  of  the  surface  shear, 
indicating  that  the  sweep  motion  is  modified. 

A  summary  plot  of  the  change  in  the  mean  surface  shear  produced  by 
three  different  widths  of  plates  at  three  different  heights  is  shown  in 
Figure  34.  Although  the  effect  of  both  plate  width  and  height  can  be 
determined  from  the  measurements,  a  large  amount  of  scatter  was  en¬ 
countered.  This  scatter  was  traced  to  a  critical  alignment  sensitivity 
of  the  plates  parallel  to  the  surface.  Figure  35  shows  tests  made  with 
the  same  plate  set  at  +2.9,  0  and  -2.9  degrees  with  respect  to  the 
surface.  As  demonstrated  by  Figure  35,  the  results  are  extremely 
sensitive  to  the  angle  of  the  plate.  Repeat  tests  for  the  plate  at  0.0 
cm  are  also  shown  in  Figure  35.  Variation  of  mean  surface  shear  stress 
as  great  as  20%  was  observed  due  to  slight  misalignment  of  the  plate 
angle  which  could  not  have  been  greater  than  ±0.5  degrees. 

The  possibility  of  surface  shear  stress  reduction  caused  by  the 
wake  behind  the  plate  was  also  investigated.  The  plate  of  3  cm  wide, 
0.4  cm  high  was  placed  with  its  trailing  edge  4  cm  upstream  of  the 
surface  sensor.  The  mean  velocity  profile  measured  above  the  surface 
sensor  at  =  8.55  m/sec  is  compared  with  result  without  the  plate  in 
Figure  36.  A  velocity  defect  was  observed  near  the  height  of  the 
plate.  The  plate  wake  showed  little  effect  on  the  flow  near  the  sur¬ 
face.  Hence,  it  can  be  concluded  that  the  wake  effect  on  the  reduction 
of  surface  shear  stress  should  be  of  secondary  importance. 

The  variations  in  the  surface  shear  stress  probability 
distribution  obtained  when  the  plate  of  3  mm  wide  and  2  mm  high  was 


placed  at  different  locations  are  compared  with  the  basic  flow  result 
in  Figure  37.  A  major  reduction  of  the  large  magnitude  fluctuations 
was  observed  when  the  plate  was  placed  at  2  cm  upstream  of  the  surface- 
hot-wire.  However,  almost  no  effect  on  the  surface  shear  stress  dis¬ 
tribution  was  seen  when  the  same  plate  was  placed  7  cm  upstream  of  the 
surface-hot-wire.  Figure  38  shows  the  surface  shear  stress  probability 
distributions  for  the  plate  of  3  mm  wide  and  8  mm  high  placed  at  var¬ 
ious  locations  upstream  of  the  surface-wire.  The  8  mm  high  plate  had 
very  little  effect  on  the  surface  shear  fluctuations. 

Based  on  the  space-time  correlation  measurements  for  U  =  S.'j.'j 
m/sec  as  shown  in  Figure  28c),  the  origin  of  the  instability  is  located 
about  0.6  cm  above  the  surface  and  6  cm  upstream  of  the  surface  wire. 
The  proposed  flow  model,  as  was  described  in  Chapter  III,  predicted 
that  the  plate  placed  on  the  passage  of  the  sweep  motion  would  have 
considerable  effect  on  the  surface  shear  stress.  This  prediction  seems 
to  be  justified  by  the  results  shown  in  Figure  37.  On  the  other  hand, 
the  results  shown  in  Figure  38  of  the  little  effect  induced  by  an  8  mm 
high  plate  implies  that  the  surface  shear  stress  is  not  affected  by 
blockage  in  the  outer  flow.  It  was  also  noted  that  the  present  tech¬ 
nique  of  blocking  the  sweep  is  effective  for  only  one  event.  Once  the 
plates  are  moved  further  upstream  their  effect  vanishes.  Thus,  stop¬ 
ping  a  single  sweep  from  reaching  the  surface  does  not  appear  to  have  a 
major  effect  on  the  overall  boundary  layer.  The  three  dimensional 
characteristics  of  the  events,  apparently  fill  in  for  the  one  blocked 
event.  Although  the  plates  are  of  major  interest  in  demonstrating  the 
general  concept  of  the  present  model  ,  they  do  not  produce  an  insight 
into  the  actual  production  of  the  sweep-burst  phenomenon. 


CHAPTER  VI 


CONCLUSIONS 

In  a  turbulent  shear  flow,  a  sudden  change  of  the  external  driving 
force  causes  a  dramatic  response  of  flow  variables  across  the  boundary 
layer.  The  response  of  the  flow  is  non-linear  due  tt-  the  viscous 
effects.  The  mixing  length  and  eddy  viscosity  models  substitute  long¬ 
time  averaged  statistical  quantities  for  the  highly  time  dependent 
variables  involved  in  a  turbulent  flow,  and  thus  mask  a  great  deal  of 
non-linear  effects.  Frandtl's  turbulent  boundary  layer  approach  of 
employing  the  averaged  Reynolds  equation  and  a  turbulence  closure  model 
has  proved  to  be  insufficient  to  predict  flow  variables  with  fluctua¬ 
tion  of  magnitude  as  great  as  the  mean  value,  such  as  the  cases  for 
surface  shear  stress  and  also  quantities  in  the  separation  region.  The 
philosophy  of  solving  first  the  time  dependent  equation  of  motion  and 
then  employing  a  stochastic  averaging  to  include  the  random  aspect  of 
the  flow  is  proposed  and  demonstrated.  This  approach  is  able  to  pro¬ 
duce  good  predictions  for  the  mean  and  time  dependent  surface  shear 
stress  under  a  turbulent  boundary  layer  over  a  smooth  plate. 

Based  on  observations  of  the  turbulent  boundary  layer  structure,  a 
physical  model  dealing  with  flow  near  the  surface  is  hypothesized. 

In  the  inner  region  of  a  turbulent  boundary  layer,  an  intermittent 
bursting  phenomenon  has  been  observed.  Also,  large  magnitude  fluctua¬ 
tions  of  the  surface  shear  stress  under  a  turbulent  boundary  have  been 
reported.  The  present  model  proposes  that  the  surface  shear  stress 
fluctuations  are  due  to  the  non-linear  response  of  the  flow  produced  by 
the  bursting  process  near  the  wall. 


It  is  postulated  that  localized  pressuie  gradients  created  in  the 
valleys  of  the  large  scale  structures  are  the  origin  of  the  bursting 
phenomenon.  This  localized  pressure  gradient  produces  the  sweep 
motion,  which  in  turn  brings  the  high-speed  fluid  from  the  outer  region 
inward  toward  the  inner  region.  The  high-speed  fluid  is  then  retarded 
near  the  wall  and  finally  ejected  outward.  The  eiectiou  of  the  low- 
speed  motion  in  turn,  triggers  the  sweep  motion  and  the  bursting  cycles 
repeat . 

A  simplified  model  was  employed  to  illustrate  the  non-linear 
response  of  surface  shear  stress  to  a  single  sweep  motion.  The  model 
treated  the  sweep  motion  as  an  impulsively  started  boundary  layer  flow 
above  a  surface.  An  approximate  time  depenent  solution  for  the  shear 
stress  at  the  wall  in  response  to  a  discrete  sweej)  velocity  was  ob¬ 
tained.  Assuming  a  Gaussian  or  a  modified  Rayleigh  probability  dis¬ 
tribution  for  the  sweep  velocity,  a  stochastic  averaging  technique  was 
developed  for  prediction  of  the  resulting  surface  shear  stress  re¬ 
sponse.  This  approach  was  shown  to  predict  the  large  magnitude  fluctu¬ 
ations  and  the  highly  skewed  probability  distribution  of  the  surface 
shear  stress. 

In  order  to  justify  the  time  dependent  model,  measurements  of  the 
probability  distribution  and  cross  rorrelat ions  of  the  turbulent  velo¬ 
city  and  surface  shear  stress  were  evaluated  toi  a  zero  pressure  gradi¬ 
ent,  turbulent  lioundary  layer.  Kxperimental  results  indicated  that  the 
probability  distributions  of  the  longitudinal  turbulent  velocity  in  the 
sublayer  region,  (0  <  ylJ^/v  <  10),  show  characteristics  similar  to  the 
probability  distribution  of  the  surface  shear  stress.  barge  magnitiicle 
fluctuations  greater  than  two  times  the  mean  sort  ace  shear  were 


observed.  Cross  ctirrc  1  nl  ion  moasiirenK-nt  s  showed  that  a  high  degree  ot 
correlation  exists  betweni  tiit  sori..ii  sliear  stress  and  the  turbulent 
velocity  in  the  inner  region.  Kroni  the  peak  i  socorrei.it  i  on  curves  near 
the  surface,  the  origin  of  the  flow  i nsLahi i i ty  can  be  identified  at  a 
non-dimensional  distance  of  approximately  (y'"  -}  IbO  above  the  surface 
and  a  non-dimensional  distance  (x-''  ~)  IbOO  upstream  of  the  surface 

sensor.  Also,  these  curves  indicate  the  appro imate  passage  of  the 
high-speed  flow  associated  with  the  sweep  motion.  It  was  also  con¬ 
cluded  that  the  convective  velocity  measured  at  the  height  near  the 
origin  of  instability  is  related  to  the  sweep  motion. 


As  implied  by  the  time  dependent  model  for  the  surface  shear 
stress,  if  the  sweep  motion  could  be  modified,  the  magnitude  of  the 
large  shear  stress  fluctuations  at  the  wall  would  be  substantially 
changed.  A  series  of  thin  (0.13  mm)  plates  of  different  widths  and 
heights  were  inserted  in  the  sublayer  to  block  the  instability  or  sweep 
motion  from  reaching  the  surface.  Reduction  of  both  the  mean  and 
fluctuating  surface  shear  was  obtained.  The  variation  in  surface  shear 
was  found  to  be  extremely  sensitive  to  slight  angle  of  attacks  of  the 
plates  ■ 
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Collecting  terms  of  like  power  of  t--'  in  Equations  (A-2),  lA-3), 
(A-4)  and  (A-5),  Equations  (19),  (20),  and  (21)  were  obtained  as 

summarized  below: 

Order  of  t"  From  Equation  Resulting  Differential  Equation 

(A-2)  2FF  ^  +  F  "  =  0 

^  o  o 

(A-2)  2^F,'  +  F,"  -  4F,  +  2F  F  ’  +  ylZGJ  ’ 

I  1  1  O  f)  1  o 
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APPENDIX  B 


DERIVATION  OF  EQUATION  132) 


To  apply  the  method  of  Rosezwejg,  a  Laplace  transform  of  f  was 


defined  as 


<t'(s,ri)  =  s  ;  f(i,t*) 


(B-l) 


Transformation  of  Equation  (30)  yields 

-  ‘  In  >n=o 


(B-2) 


where  s  is  very  large.  For  small  s,  transformation  of  Equation  (31) 


becomes 


_  f  ~  _f  >  foo') 

^  3n  ^n=o  -  ^  ^  ^ 


(B-3) 


Following  Rosenzweig,  an  approximate  expression  will  be  assumed  of 


the  form 

f  ?i  ^  . 

an  ’n=o  ' 

For  small  s  this  reduces  to 


(B-4) 


.  ( 91 )  =  :ii2 , 

(an  Vo  ••• 


while  for  large  s 


.  (  9$  )  .  ,  3I/2  .  J  .  ,-1/2  ^ 

9n  -^0=0  1  1  (  2  2  ■ 


(B-3) 


(B-6) 


where  a^,  and  a^  were  determined  such  that  Equation  (B-5)  agree 
with  Equation  (B-3),  and  Equation  (B-6)  agree  with  Equation  (B-2).  We 


•1  =  -^/2 


a^  =  2f'^(») 
32  =  f'^(®) 


have 


Thus  for  any  the  approximate  form  for  ^(p-Q  obtained  by 

^J\  =  -r'^  ( .  •  ) 

9niri=o  '  ^  ^ 


a  (s+a  ) 
=  -L  ‘  {  -i-  -  -  ^- 


iVs+a^ 


=  -L  •  {  — ^ 


a  a 

+  -  } 


Vs+a^  sVs+a^ 

=  -1.-1  (  _!i_  }  -  L‘*  {  — ^  } 

Vs+a^  s^/s+a^ 

=  exp  {-2f erf  {^f2t*f'(“)} 

4nt^' 
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APPENDIX  D 

EVALUATION  OF  THE  MAXIMUM  SPACE-TIME  CORRELATIONS 


Data  for  cross-correlation  measurements  between  a  boundary-layer- 

hot-wire  placed  at  coordinate  (x,y)  and  a  surface-wire  located  at 

coordinate  (0.0, 0.0)  are  presented.  Positive  x  is  designated  as 

upstream  of  the  surface-wire.  (R  u)  is  the  maximum  cross-correla- 

Tw  max 

tion  value  obtained  at  a  maximum  time  delay,  At  ,  between  the 

^  max 

boundary-layer-  and  surface-hot-wire  signals.  The  boundary-layer-hoL- 
wire  signal  was  delayed  with  respect  to  the  surface  wire. 


a) 

10.20  m/sec,  R„  =  4050 
o 

b) 

It 

9.36  m/sec,  Kq  =  3500 

c) 

8.55  m/sec,  Rg  =  3080 
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0. 

258 

- 

-4, 

.35 

0, 

.  127 

0. 

,309 

-4 

.50 

-5. 

75 

1 . 

219 

0. 

252 

) 

10 

-4 

.35 

0 

.203 

0 

.326 

-4 

.50 

-4 

.35 

0 

.330 

0, 

.356 

-3 

.  75 

-6. 

45 

0. 

127 

9  . 

114 

-9 

00 

-4 

.  35 

0 

.457 

0 

.366 

-3 

.50 

-6 . 

45 

0. 

254 

0. 

304 

-8 

.56' 

-4 

.35 

0 

.584 

0, 

.326 

-2 

.25 

-6, 

45 

0. 

381 

0. 

335 

-  7 

50 

-4 

.35 

0 

.711 

0 

.414 

-0 

.60 

-6. 

45 

0, 

508 

0, 

,149 

-  7 

''0 

.35 

0 

.838 

0, 

.296 

1 

.00 

-6. 

45 

0, 

635 

0, 

317 

-0 

03 

-4 

.35 

0 

.  965 

0 

.263 

1 

.00 

-6. 

45 

0, 

762 

0. 

448 

-4 

00 

-4 

.35 

1 

.092 

0 

.270 

1 

.  50 

-6. 

45 

0, 

889 

0. 

263 

-4 

00 

-4 

.  35 

1 

.219 

0, 

.  20-1 

2 

.40 

-6. 

45 

i  . 

016 

0. 

2  7  4 

-3 

0(' 

-6. 

45 

1 , 

14.3 

0  , 

245 

- 

-6. 

45 

1 , 

270 

0. 

2.30 

-0 

50 

-6 . 

95 

0 

0  76 

0, 

314 

-10 

50 

-t). 

95 

0, 

•)Q7 

0. 

327 

-8 

50 

-6 

95 

0 

330 

0  , 

409 

“  7 

75 

-6, 

95 

0 

457 

0 

305 

-7 

00 

-6 

95 

0 

584 

265 

-  4 

00 

-6 

95 

0 

711 

0 

349 

-k 

00 

-6 

95 

0 

838 

0 

2  79 

-4 

00 

9  5 

0 

965 

0 

2  i  2 

-  4 

00 

-(> 

93 

1 

09.’ 

0 

280 

-1 

(!l) 

95 

1 

219 

0 

2  7  1 

-0 

Vt 

76 


rj  U  ^  8.5'j  in/ set,  K,,  -  i()80 

OC 


X  (cm)  y  (cm) 


(H.  ii) 
iW  max 


At 


max 
lO'  '  set 


X  (cm)  y  (cm) 


iK  u) 

Tw  max 


At 


max 


10  ^  sec 


10.8 

0. 

127 

0. 

168 

24. 

50 

8. 

3 

0. 

051 

0. 

175 

16. 

00 

10.8 

0. 

254 

0. 

222 

21 

75 

8. 

3 

0. 

102 

0. 

163 

16. 

50 

10.8 

0. 

381 

0. 

227 

25 

00 

8. 

3 

0. 

229 

0 

218 

21 

25 

10.8 

0. 

445 

0. 

200 

24. 

75 

8. 

3 

0. 

356 

0. 

290 

21  . 

50 

10.8 

0. 

508 

0. 

261 

24. 

00 

8 

3 

0. 

419 

0 

305 

21 

50 

10.8 

0. 

584 

0. 

239 

24. 

50 

8. 

3 

0. 

483 

0. 

268 

21 

50 

10.8 

0 

686 

0. 

296 

28 

00 

8 

3 

0 

559 

0 

260 

21 

50 

10.8 

0. 

762 

0. 

284 

26 

00 

8 

3 

0 

660 

0 

245 

21 

50 

10.8 

0 

889 

0. 

270 

25 

25 

8 

3 

0 

737 

0 

283 

21 

50 

10.8 

1. 

016 

0. 

292 

25 

00 

8 

3 

0 

864 

0 

270 

20 

50 

8 

3 

0 

991 

0 

257 

21 

75 

10.3 

0 

102 

0 

170 

21 

50 

8 

3 

1 

245 

0 

206 

22 

50 

10.3 

0 

229 

0 

196 

22 

00 

10.3 

0 

356 

0 

188 

24 

00 

7 

6 

0 

076 

0 

161 

15 

65 

10.3 

0 

419 

0 

188 

22 

25 

7 

6 

0 

178 

0 

172 

17 

65 

10.3 

0 

483 

0 

200 

23 

25 

7 

6 

0 

254 

0 

239 

17 

00 

10.3 

0 

559 

0 

208 

23 

25 

7 

6 

0 

381 

0 

244 

18 

20 

10.3 

0 

660 

0 

310 

25 

00 

7 

6 

0 

508 

0 

278 

19 

00 

10.3 

0 

737 

0 

240 

23 

50 

7 

6 

0 

635 

0 

269 

20 

25 

10.3 

0 

864 

0 

220 

23 

.50 

7 

.6 

0 

762 

0 

318 

21 

00 

10.3 

0 

991 

0 

198 

26 

75 

7 

6 

0 

889 

0 

287 

17 

00 

7 

6 

1 

016 

0 

242 

18 

50 

8.8 

0 

076 

0 

179 

18 

50 

8.8 

0 

127 

0 

153 

22 

00 

7 

0 

0 

051 

0 

208 

14 

20 

8.8 

0 

254 

0 

218 

22 

00 

7 

.0 

0 

152 

0 

211 

17 

15 

8.8 

0 

.381 

0 

286 

22 

.00 

7 

.0 

0 

229 

0 

236 

16 

60 

8.8 

0 

445 

0 

287 

21 

50 

7 

0 

0 

356 

0 

267 

18 

00 

8.8 

0 

.508 

0 

288 

21 

.50 

7 

.0 

0 

483 

0 

.286 

14 

50 

8.8 

0 

.584 

0 

454 

22 

00 

7 

.0 

0 

610 

0 

290 

15 

50 

8.8 

0 

.686 

0 

.264 

21 

.50 

7 

.0 

0 

.762 

0 

.250 

15 

.00 

8.8 

0 

.762 

0 

.238 

22 

.00 

7 

.0 

0 

864 

0 

191 

16 

.75 

8.8 

0 

.889 

0 

254 

22 

.50 

7 

.0 

0 

,991 

0 

.241 

16 

.00 

8.8 

1 

.016 

0 

.271 

21 

.75 

7 

.0 

1 

.245 

0 

.176 

16 

.25 

8.8 

1 

.270 

0 

.200 

22 

.00 

■ar 


r  ■ 


c)  =  8.55  m/sec,  =  3080  (continued) 


X  (cm) 

y  (cm) 

(R,  u) 

Tw  max 

At 

max 

X  (cm) 

y  (cm) 

(R  u ) 

TW  max 

At 

max 

10  ^  sec 

10  ^  .s  e  c 

6 . 4 

0.076 

0.231 

12.60 

4.7 

0.076 

0.233 

10.15 

6.4 

0.178 

0 . 205 

14.80 

4.7 

0.  178 

0.290 

12.50 

6.4 

0.254 

0.218 

15.40 

4.7 

0.254 

0 . 322 

11.50 

6 . 4 

0.381 

0 . 305 

14.00 

4.7 

0.381 

0.320 

12.00 

6 . 4 

0.508 

0.320 

15.62 

4.7 

0.508 

0 . 320 

11.20 

6.4 

0.635 

0.290 

15.70 

4.7 

0.635 

0.290 

12.10 

6.4 

0.762 

0.260 

16.00 

4.7 

0,762 

0.260 

13.25 

6.4 

0.889 

0.283 

16.00 

4.7 

0.889 

0.287 

13.00 

6.4 

1.016 

0.218 

18.40 

4.7 

1.016 

0.240 

14.45 

6.4 

1.270 

0.196 

19.20 

4.7 

1  .  270 

0.190 

15.00 

5.8 

0.051 

0.239 

12.45 

4.1 

0 . 064 

0 . 250 

10.00 

5.8 

0.165 

0.239 

12.50 

4.1 

0.  165 

0.275 

10.65 

5.8 

0.241 

0.240 

13.75 

4.1 

0.241 

0.297 

11.80 

5.8 

0.368 

0.297 

15.40 

4.1 

0.368 

0.314 

11.10 

5.8 

0.495 

0.281 

16.00 

4.1 

0.495 

0.326 

12.45 

5.8 

0.622 

0.272 

15.20 

4.  1 

0.622 

0.293 

11.90 

5.8 

0.749 

0.225 

15.30 

4. 1 

0.749 

0.269 

12.08 

5.8 

0.876 

0.226 

16.85 

4.1 

0.876 

0.2  70 

11.20 

5.8 

1  .003 

0.246 

17.65 

4.1 

1.003 

0.205 

13.55 

5.8 

1.257 

0.177 

18.00 

4. 1 

1.257 

0.  too 

13.6(1 

5.2 

0.076 

0.200 

12.50 

3.4 

0.051 

0.309 

9.30 

5  2 

0.254 

0.251 

13.75 

3.4 

0.127 

0.325 

9.17 

5.2 

0.381 

0.293 

16.00 

3.4 

0.254 

0.376 

10,16 

5.2 

0.508 

0.292 

15.00 

3.4 

0.381 

0.494 

9.53 

5.2 

0.635 

0.258 

14.25 

3.4 

0.508 

0.266 

9.84 

5.2 

0.762 

0.286 

16.00 

3.4 

0.762 

0.299 

10.50 

5.2 

0.889 

0.248 

17.00 

3.4 

1  .270 

0.189 

11  .73 

5.2 

1.016 

0 . 249 

16.00 

78 


c)  III/ ::i'(  ,  (ronl  i  iiui'il  J 


X  (cm) 

y  (cm) 

(R  u ) 

Tw  max 

max 

lo’^  si-c 

X  (  cm  ) 

y  (cm) 

(i<  u) 

Tw  max 

At 

max 

lO"^  se( 

2.4 

0 . 0,3 1 

0 . 345 

6 . 46 

0.15 

0.025 

0.685 

0.75 

2.4 

0.  127 

0.424 

7.64 

0.15 

0.076 

- 

1.025 

2.4 

0.2.34 

0.370 

8.40 

0.15 

0.127 

- 

2.  kb 

2.4 

0.381 

0.411 

8.62 

0.15 

0.178 

0.550 

3.00 

2.4 

0.508 

0 . 350 

8.40 

0.15 

0.279 

0.449 

3.40 

2.4 

0.762 

0.330 

9.60 

0.15 

0.356 

0.378 

4.00 

2.4 

1.270 

0.209 

11.80 

0.15 

0.508 

0.322 

5.00 

0.15 

0.635 

0.333 

5.50 

1.7 

0.127 

0.496 

5.50 

0.15 

0.762 

0.292 

5.75 

1.7 

0.254 

0.400 

6.50 

0. 15 

0.889 

0.272 

6.50 

1.7 

0.381 

0.340 

7.00 

0.15 

1.016 

0.246 

7.50 

1.7 

0.508 

0.342 

7.50 

1.7 

0.635 

0.301 

8.75 

-0.15 

0.051 

0.795 

0.65 

1.7 

0.762 

0.292 

9.00 

-0.15 

0.  102 

0.693 

1.00 

1.7 

0.889 

0.290 

9.50 

-0.15 

0.152 

0.639 

1  ..30 

1.7 

1.016 

0.253 

9.50 

-0.15 

0.203 

0.517 

2.25 

-0.15 

0.254 

0.452 

2.50 

1.3 

0.051 

0.462 

.3.74 

-0.15 

0.305 

0.446 

3.25 

1.3 

0.254 

0.441 

4.92 

-0.15 

0.381 

0.390 

3.60 

1.3 

0.381 

0.393 

4.72 

-0.15 

0.508 

0.382 

5.20 

1.3 

0.508 

0.332 

5.87 

-0.15 

0.762 

0.278 

6.15 

1.3 

0.762 

0.325 

8.19 

-0.15 

1.270 

0.294 

10.50 

0.7 

0.51 

0.590 

1.65 

-0.75 

0.025 

0.774 

-0.55 

0.7 

0. 102 

0.740 

1.85 

-0.75 

0.076 

0.716 

0.40 

0.7 

0.203 

0.457 

3.10 

-0.75 

0. 127 

0.617 

1.30 

0.7 

0.305 

0.409 

4. 15 

-0.75 

0.178 

0.505 

0.7 

0.381 

0.432 

3.90 

-0.75 

0.229 

0.465 

1.50 

0.7 

0.508 

0.436 

5.70 

-0.75 

0.279 

0.448 

2.50 

-0 . 75 

0 . 356 

0.388 

2.99 

0.9 

0.635 

0.309 

6.00 

-0.75 

0.483 

0.352 

3.75 

0.9 

0.762 

0.335 

7.00 

-0.75 

0.635 

0.338 

6.50 

0.9 

0 . 889 

0.249 

9.00 

-0.75 

0.762 

0.310 

6.00 

0.9 

0.143 

0.240 

10.00 

-0.75 

0.889 

0.313 

6.25 

-0.75 

1.016 

0.288 

7.00 

79 


c)  =  8.55  in/sec,  =  3080  (continued) 


X  (cm)  y  (cm) 

Tw  max  max 


X  (cm)  y  (cm)  (R  u)  At 

TW  max  max 


10-3 

sec 

lO"^ 

St’i 

-2. 

,1 

0. 

,025 

0, 

,548 

-3. 

,90 

-3. 

,85 

0. 

,  102 

0. 

250 

-8. 

,50 

-2. 

,  1 

0, 

,076 

0. 

.444 

-3. 

,05 

-3. 

,85 

0. 

178 

0. 

263 

-6. 

,25 

-2. 

.  1 

0. 

.127 

0. 

.533 

-1 . 

.90 

-3. 

.85 

0. 

,254 

0. 

,316 

-4. 

,00 

-2. 

,1 

0. 

.178 

0, 

.432 

-1 . 

,80 

-3. 

.85 

0. 

.381 

0. 

269 

-3. 

,25 

-2. 

.  1 

0, 

,229 

0. 

.488 

-0. 

25 

-3, 

,85 

0. 

,508 

0. 

315 

-2. 

50 

-2, 

.1 

0, 

.356 

0. 

.384 

-0. 

.25 

-3. 

.85 

0. 

,635 

0. 

356 

-4. 

00 

-2. 

,  1 

0, 

,483 

0. 

.370 

0. 

,0 

-3. 

,85 

0. 

762 

0. 

.300 

-1 . 

25 

-3. 

.85 

0. 

,889 

0. 

.291 

-0. 

25 

-2. 

.25 

0, 

,635 

0. 

.295 

2. 

.25 

-3. 

.85 

1  , 

,016 

0. 

302 

1 , 

00 

-2. 

.25 

0. 

,762 

0. 

.279 

2. 

.50 

-3. 

.85 

1  , 

.270 

0. 

.224 

I . 

.50 

-2. 

,25 

0. 

.889 

0. 

.286 

3. 

.00 

-2. 

,25 

1 . 

,016 

0. 

.261 

3. 

.25 

-4. 

,35 

0. 

.051 

0. 

,252 

-8. 

.25 

-2. 

25 

1 . 

270 

0. 

,221 

7. 

.00 

-4, 

.35 

0. 

127 

0, 

.326 

-6 

.00 

-4. 

.35 

0. 

.203 

0, 

,300 

-7 

.50 

-2. 

6 

0. 

,051 

0. 

.519 

-5. 

,00 

-4, 

.35 

0 

.  330 

0, 

.331 

-4 

.00 

-2. 

6 

0. 

,  102 

0, 

,434 

-4. 

.00 

-4, 

.35 

0, 

.457 

0, 

.344 

-3 

.50 

-2, 

,6 

0. 

.152 

0. 

,474 

-3. 

.20 

-4. 

.35 

0. 

.584 

0, 

.332 

-1 

.■;)0 

-2. 

,6 

0. 

,203 

0. 

.477 

-3, 

.50 

-4, 

.35 

0. 

.711 

0, 

,386 

-0 

.50 

-2. 

.6 

0. 

.254 

0, 

,421 

-2, 

.70 

-4, 

.35 

0, 

.838 

0, 

,241 

0 

.50 

-2. 

,6 

0. 

.381 

0. 

,406 

-1 , 

.50 

-4, 

.35 

0, 

.965 

0. 

.313 

-0. 

.25 

-2. 

,6 

0. 

,508 

0. 

,366 

-0, 

,90 

-4, 

.35 

1 , 

,219 

0. 

.237 

1 , 

.00 

-2. 

,75 

0. 

,711 

0, 

.285 

1  , 

,50 

-5. 

.25 

0, 

,  102 

0. 

,291 

-8, 

50 

-2. 

,75 

0. 

,838 

0 

.340 

1 . 

.20 

-5. 

.25 

0. 

,254 

0. 

288 

-6. 

00 

-2. 

,75 

0. 

,965 

0 

.242 

3. 

,00 

-5. 

.25 

0. 

,381 

0. 

,288 

-4, 

,50 

-2. 

.75 

1 , 

.219 

0, 

.250 

6, 

.25 

-5. 

.25 

0. 

.508 

0. 

300 

-4. 

,00 

-5. 

.25 

0. 

.635 

0. 

,332 

-4. 

,50 

-3. 

.2 

0, 

.025 

0 

.328 

-6, 

.65 

-5, 

.25 

0, 

.762 

0. 

329 

-2, 

00 

-3. 

.2 

0 

.076 

0 

.521 

-5, 

.85 

-5, 

.25 

0, 

.889 

0. 

280 

-3. 

,25 

-3 

.2 

0 

.127 

0 

.516 

-4 

.  10 

-5. 

.25 

1 , 

,016 

0  , 

274 

-  4 . 

,50 

-3. 

.2 

0 

.  178 

0 

.523 

-3, 

.80 

-5, 

.25 

1 . 

,  270 

0. 

,209 

1 . 

,50 

-3, 

0, 

.229 

0 

.486 

-3 

,00 

-3. 

,2 

0, 

.356 

0 

.411 

-1 

.40 

-3 

.2 

0 

cc 

0 

.420 

-0 

.50 

c)  11  ::  8 

a. 

.55  iit/sec,  =  3080  ( c onl  iniicd  ) 

X  (cm) 

y  (tm) 

(R  u) 

At 

Tw  max 

max 

-3 

10  '  SIM 

-5.75 

0 . 05  1 

0.2  79 

1  .05 

-5.75 

0.203 

0.297 

-7.25 

-5.75 

0.330 

0.357 

-6.25 

-5.75 

0.457 

0 .  i05 

-4.50 

82 


Table  2.  Evaluation  of  Equation  (32) 


Nondimensional  Time 


Nondimensional  Shear  Stress 

-I 

3r|ln=0 


0.0036 

0.0082 

0.0145 

0.0227 

0.0327 

0.0510 

0.0655 

0.0907 

0.1098 

0.1533 

0.2041 

0.2622 

0.3450 

0.4391 

0.5671 

0.7896 

0.9291 

1.0800 

1.2421 

1.6005 

2.0042 

2.5008 

3.2121 

3.8334 

4.7035 

5.7344 

6.9477 

8.1026 

9.0732 


13.244 

8.829 

6.622 

5.298 

4.415 

3.532 

3.117 

2.649 

2.409 

2.039 

1.768 

1.562 

1.363 

1.211 

1.069 

0.914 

0.8475 

0.7919 

0.7449 

0.6707 

0.6162 

0.5724 

0.5342 

0.5137 

0.4962 

0.4847 

0.4774 

0.4740 

0.4722 


83 


Table  3.  II  i  st  ogr.ini  and  esliiiiated  probability  'len;;ity  (Or  f.  Values 
'  between  * ~  ■  . . 


0  anil  t" 


4.0  were  used. 


n'  aI- 

o  o 

N  ^00  ,  Af’  0.1 
»>  o 


window 

N  . 

1 

P 

1 

0.5 

. 

0.6 

95.0 

4.75 

0.6 

- 

0.7 

32.0 

1.60 

0.7 

- 

0.8 

20.0 

1.00 

0.8 

- 

0.9 

12.5 

0.625 

0.9 

- 

1.0 

7.8 

0.390 

1.0 

- 

1.1 

5.9 

0.295 

1.1 

- 

1.2 

4.8 

0.240 

1.2 

- 

1.3 

3.4 

0.170 

1.3 

- 

1.4 

2.7 

0.135 

1.4 

- 

1.5 

2.0 

0.100 

1.5 

- 

1.6 

1.9 

0.095 

1.6 

- 

1.7 

1.4 

0.070 

1.7 

- 

1.8 

1.0 

0.050 

1.8 

- 

1.9 

0.95 

0.0475 

1.9 

- 

2.0 

0.85 

0.0425 

2.0 

- 

2. 1 

0.74 

0.0370 

2.1 

- 

2.2 

0.66 

0.0330 

2.2 

- 

2.3 

0.64 

0.0320 

2.3 

- 

2.4 

0.56 

0.0280 

2.4 

- 

2.5 

0.47 

0.0235 

2.5 

- 

2.6 

0.38 

0.0190 

2.6 

- 

2.7 

0.34 

0.0170 

2.7 

- 

2.8 

0.30 

0.0150 

2.8 

- 

2.9 

0.28 

0.0140 

2.9 

- 

3.0 

0.26 

0.0130 

'  ) 

O  ! 

w i ndow 

N. 

I 

P. 

1 

3.0 

- 

3.1 

0.25 

0.0125 

3.1 

- 

3.2 

0.23 

0.0115 

3.2 

- 

3.3 

0.20 

0.0100 

3.3 

- 

3.4 

0.20 

0.0100 

3.4 

- 

3.5 

0.18 

0.0090 

3.5 

- 

3.6 

0.15 

0.0075 

3.6 

- 

,3.7 

0.13 

0 . 0065 

3.7 

- 

3.8 

0.1 1 

0.0055 

3.8 

- 

3.9 

0.  10 

0.0050 

3.9 

- 

4.0 

0.09 

0.0045 

4.0 

- 

4. 1 

0.06 

0.0030 

4.1 

- 

4.2 

0.04 

0.0020 

4.2 

- 

4.3 

■,v 

Note:  Too  liLt.lr  to  be  counted 


Table  4.  Parameters  used  in  the  stochastic  averaging  computation  and  the  mean  surface  shear  stress 
obtained. 
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Table  5.  Static 
model . 

pres.sure  distribution 
X  wa.s  measured  from 

along  the  ceritc 
the  entrance  of 

rline  of  the  test 
the  test  section. 

Freestream 

Velocity 

U  10.7 

00 

m/sec  U 

00 

=  9.76  m/sec 

U  =  8.55  m/sec 

x-dibtance 

Ap 

Ap 

Ap 

(cm) 

(tnmHg) 

(mmHg) 

(mmHg) 

120.40 

+0.011 

+0.008 

+0.006 

130.65 

-0.001 

-0.0016 

-0.002 

135.72 

+0.002 

+0.001 

+0.000 

140.77 

+0.000 

-0.000 

-0.001 

145.85 

+0.000 

-0.000 

-0.001 

151.01 

+0.002 

+0.0006 

-0.000 

156.01 

+0.003 

+0.001 

+0.000 

161.09 

+0.001 

+0.000 

-0.000 

166. 12 

+0.002 

+0.001 

+0.000 

171.20 

+0.001 

-0.001 

-0.000 

176.28 

+0.002 

+0.001 

+0.000 

181.36 

+0.003 

+0.002 

+0.0009 

186.44 

+0.002 

+0.001 

+0.000 

191.52 

+0.002 

+0.001 

+0.000 

196.60 

-0.001 

-0.002 

-0.002 

201.68 

-0.003 

-0.0026 

-0.0024 

206.76 

-0.007 

-0.006 

-0.006 

209.30 

+0.003 

+0.0015 

+0.0005 

211.82 

-0.016 

-0.014 

-0.012 

214.36 

-0.025 

-0.021 

-0.018 

216.90 

-0.031 

-0.026 

-0.021 

218.21 

-0.041 

-0.036 

-0.028 

219.44 

-0.055 

-0.046 

-0.038 

220.07 

-0,064 

-0.055 

-0.044 

220.71 

-0.072 

-0.060 

-0.050 

220.99 

-0.077 

-0.065 

-0.052 

221.34 

-0.084 

-0.070 

-0.058 

222.80 

-0.082 

-0.071 

-0.056 

223.415 

-0.069 

-0.059 

-0.050 

224.705 

-0.065 

-0.054 

-0.045 

227.245 

-0.050 

-0.043 

-0.035 

229.785 

-0.037 

-0.031 

-0.026 

232.325 

-0.028 

-0.024 

-0.020 

234.883 

-0.020 

-0.016 

-0.014 

237.405 

-0.013 

-0.011 

-0.010 

239.945 

-0.003 

-0.004 

-0.003 

242.485 

+0.007 

+0.006 

+0.004 

245.025 

+0.024 

+0.020 

+0.015 

247 . 565 

+0.042 

+0.035 

+0.028 

250. 105 

+0.058 

+0.048 

+0.039 

252.66 

+0.067 

+0.055 

+0.041 

255.20 

to. 076 

* 

+  0.0(>4 

+0.050 

. . 
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Table  5.  (continued) 


Freestream 

Velocity 

U  =  10.7  m/sec 

=  9.76  m/sec 

U  =  8.55  r 

CO 

x-distance 

Ap 

Ap 

Ap 

(cm) 

(mmHg) 

(mmHg) 

(mmHg) 

257.74 

+0.083 

+0 . 069 

+0.054 

260.30 

+0.093 

+0.077 

+0.061 

265.43 

+0.097 

+0.081 

+0 . 064 

267.95 

+0.100 

+0.085 

+0.066 

270.46 

+0.104 

+0.087 

+0.068 

272.98 

+0.107 

+0.089 

+0.071 

275.52 

+0.111 

+0.091 

+0.072 

278.11 

+0.112 

+0.094 

+0.074 

280.67 

+0.114 

+0.096 

+0.074 

283.19 

+0.116 

+0.099 

+0.077 
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Table  6.  Boundary  layer  quantities. 


Uoo 

(m/sec) 

6 

(cm) 

6* 

(cm) 

0 

(cm) 

H 

10-3 

T 

w 

(N/m^) 

U 

T 

(m/ sec) 

10.20 

6.350 

0.991 

0.718 

1.380 

4050 

3.225 

0.1713 

0.4147 

9.36 

6.477 

0.932 

0.690 

1.350 

3500 

3.350 

0.1439 

0.3800 

8.55 

6.528 

0.879 

0.664 

1.324 

3080 

3.450 

0.1243 

0.3532 

Table  7.  Measured  mean  velocity  distribution  over  the  flat  plate. 


a)  =  10.20  m/sec 


y  y/6  yu^/v  u  u/u^  u/u^ 


cm 

m/sec 

0.018 

0.0030 

4.00 

1.73 

0.170 

4.17 

0.023 

0.0036 

5.19 

2.27 

0.225 

5.47 

0.028 

0.0044 

6.31 

2.48 

0.243 

5.98 

0.031 

0.0049 

6.99 

3.08 

0.302 

7.42 

0.045 

0.0071 

10.14 

3.59 

0.352 

8.65 

0.058 

0.0090 

13.08 

4.23 

0.415 

10.20 

0.083 

0.0131 

18.71 

4.65 

0.451 

11.21 

0.100 

0.0157 

22.54 

5.04 

0.494 

12.15 

0.134 

0.0211 

30.21 

5.25 

0.515 

12.66 

0.160 

0.0252 

36.08 

5.42 

0.531 

13.07 

0.210 

0.0331 

47.35 

5.76 

0.565 

13.89 

0.262 

0.0413 

59.08 

5.79 

0.568 

13.96 

0.388 

0.0611 

87.49 

6.28 

0.616 

15.14 

0.516 

0.0811 

116.35 

6.55 

0.642 

15.79 

0.642 

0.1011 

144.77 

6.79 

0 . 666 

16.37 

0.770 

0.1213 

173.63 

6.92 

0.678 

16.68 

1.024 

0.1613 

230.90 

7.10 

0.696 

17.12 

1.278 

0.2013 

288.18 

7.34 

0.  720 

17.70 

2.040 

0.3213 

460.00 

8.10 

0.794 

19.53 

2.548 

0.4013 

574.57 

8.36 

0.820 

20.17 

3.310 

0.5213 

746.40 

8.65 

0.848 

20.86 

3.818 

0.6013 

861.00 

9.15 

0.897 

22.06 

4.580 

0.7213 

1032.80 

9.44 

0.926 

22.78 

5.088 

0.8013 

1147.33 

9.69 

0.950 

23.37 

5,586 

0.8800 

1259.00 

9.82 

0.963 

23.69 

5.842 

0.9200 

1317.00 

10.00 

0.980 

24,11 

6.092 

0.9594 

1373.36 

10.09 

0.989 

24.33 

6.350 

1.0000 

1431.00 

10.20 

1.000 

24.60 

8<J 


Ti.hli-  7. 

(  C'UIlL  1  lll.Cil  ) 

= 

JO 

9,36  iii/sfH- 

y 

y/6 

vlJ  /v 

"  T 

U 

u/u^ 

cm 

m/ser 

0.022 

0.0034 

4.54 

1  .75 

0. 188 

4.60 

0.027 

0.0042 

5.57 

2.00 

0.215 

5.28 

0.033 

0.0051 

6.82 

2.28 

0 . 245 

6.02 

0.037 

0.0057 

7.64 

2.91 

0.312 

7.66 

0.047 

0.0073 

9.71 

3.32 

0.357 

8.74 

0.063 

0.0097 

1 3 . 02 

3 . 66 

0.393 

9.63 

0.073 

0.0113 

15.08 

3.91 

0.420 

10.29 

0.087 

0.0134 

17.97 

4, 10 

0.440 

10.80 

0.111 

0.0171 

22.94 

4.52 

0.485 

11.39 

0.163 

0.0252 

36 . 68 

4.95 

0.531 

13.03 

0.213 

0.0329 

44.01 

5.24 

0.562 

13.78 

0.263 

0.0406 

54.34 

5.43 

0.580 

14.28 

0.391 

0.0604 

80.79 

5.65 

0.606 

14.87 

0.516 

0.0797 

106.62 

5.79 

0.622 

15.25 

0.642 

0.0991 

132.66 

5.96 

0.639 

16.68 

0.769 

0.1187 

158.90 

6.20 

0 . 665 

16.31 

1.023 

0.1579 

211.38 

6 . 36 

0.680 

16.74 

1.530 

0.2362 

316.15 

6.89 

0.740 

18.14 

2.041 

0.3154 

422.15 

7.25 

0.778 

19.09 

2.546 

0.3931 

526.08 

7.51 

0.803 

19.77 

3.053 

0.4714 

630.84 

7.89 

0.843 

20.76 

3.814 

0.5889 

788.10 

8.38 

0.896 

22.06 

4.575 

0.7063 

945 . 33 

8.72 

0.936 

22.95 

5.084 

0.7849 

1050.51 

8.80 

0.940 

23.14 

5.690 

0.8785 

1175.73 

9.03 

0.965 

23.76 

6.092 

0.9450 

1258.80 

9.20 

0,983 

24.20 

6.477 

1.0000 

1337.52 

9 . 36 

1  .000 

24.62 
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Table  8,  ITirbii  leiice  i  ill  i-iit,  i  I  y  ineasii  t  emeu  t  s  . 


a)  -  10.  ZO  III/ sec 


V 

y/6 

Turbulent  Intensity 

(cnO 

(%) 

0 . 008 

0 . 00 1 

29 . 0 

0.020 

U.OO,} 

28.2 

0 . 038 

0.00b 

25.3 

0.051 

0.008 

23.2 

0.076 

0.012 

20.5 

0.127 

0.020 

16.  / 

0.190 

0.030 

14.3 

0.254 

0 . 040 

13.4 

0.381 

0.060 

11.2 

0.508 

0.080 

10.3 

1.016 

0.160 

9.0 

1.524 

0.240 

8.0 

2.032 

0.320 

7.6 

2.540 

0.400 

6.2 

3.302 

0.520 

5.1 

4.318 

0.680 

3.7 

5.080 

0.800 

3.0 

5.842 

0.920 

2.0 

a)  Ug,  =  10.20  m/sec 


y 

(cm) 

y/6 

Turbulent  Intensity 
(%) 

0.008 

0.001 

27.7 

0.020 

0.003 

29.0 

0.038 

0.006 

24.6 

0.051 

0.009 

20.8 

0.076 

0.01 1 

20.1 

0.127 

0.019 

17.3 

0.190 

0.029 

14.8 

0.254 

0.039 

13.3 

0.381 

0.058 

11.8 

0.508 

0.078 

10.7 

1.016 

0.157 

9.0 

1.524 

0.235 

7.9 

2.032 

0.313 

7.3 

2.540 

0.392 

6.2 

3.302 

0.509 

5.0 

4.318 

0  666 

3.9 

5 . 080 

0.784 

3.0 

5.715 

0.882 

1.7 

0.008 

0.001 

27.3 

0.020 

0.003 

32.0 

0.038 

0.005 

24.6 

0.051 

0.007 

20.4 

0.076 

0.011 

20.3 

0.127 

0.019 

17.1 

0. 190 

0.029 

15.2 

0.254 

0.039 

13.7 

0.381 

0.058 

11.6 

0.508 

0.077 

11.1 

1.016 

0. 155 

9.  1 

1.524 

0.233 

8.2 

2.032 

0.311 

7.3 

2.540 

0.389 

6.5 

3 . 302 

0.505 

5.1 

4.318 

0.661 

4.3 

5.082 

0.778 

2.4 

5.715 

0.875 

1.6 
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Table  9.  Intermittency  measurements. 


U  =  10. 
00 

,20  m/sec 

00 

.36  m/sec 

-  8.55 

m/ser 

y/6 

Y 

y/6 

Y 

y/6 

Y 

0.400 

0.97 

0.387 

0.98 

0.395 

0.98 

0.500 

0.96 

0.492 

- 

0.494 

0.92 

0.600 

0.95 

0.589 

0.93 

0.585 

0.87 

0.700 

0.93 

0.686 

0.92 

0.684 

0.83 

0.800 

0.77 

0.782 

0.51 

0.782 

0.77 

0.900 

0.18 

0.879 

0.23 

0.881 

- 

1.003 

0.03 

0.976 

0.14 

0.980 

0.08 

Probability ,  P(  r 


Figure  I4b.  Predicted  probability  distribution;  Gaussian  distribution 
Uc=4.3  m/sec,  cr=1.0  ni/sec. 
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Figur 


ability,  P(  r*) 


Probability  ,  P  (  t 


Figure  Predicted  probability  distribution;  Gaussian  distribution 

Uc=5.1  m/sec,  0=1.5  m/sec. 


ProbobHity  ,  P(t 


Figure  I4g.  Predicted  probability  distribution;  Rayleigh  distribution 
(Uc)min  =  3.2  m/sec,  c  =  1.0  m/sec. 


•  J 

■J 


Probability,  P( 
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Figure  lAh.  Predicted  probability  distribution;  Rayleigh  distribution 
(Uc)min  =  3.0  m/sec,  c  =  1.3  m/sec. 
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Figure  I4i.  Predicted  probability  distribution;  Rayleigh  distribution 
(Uc)inin  =  2.5  m/sec,  c  =  1.8  m/sec. 
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Surfoce  sheer  sfess  ,  ,  N/m‘ 


f’redi  i:tc<l  probability  <li  st  r  i  but  i  on  ;  Rdyloi>>h  d  i  si  r  i  but  i  on 
(Ur)min  -  :i.O  m/sei-  t  =  1.9  tn/sor. 
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Figure  17.  Static  pressure  distribution  along  the  model  surface. 


cm 


Figure  18.  Mean  velocity  profiles  measured  for  several  lateral 

locations  along  the  test  model.  =  12.80  m/sec  was 
measured  nt  x  =  0.0  cm. 
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F i gu  re  19. 


One  parameter  ecpialion  of  Bell  (1979)  tor  the  skin 
friction  evaluation  over  a  zero  pressure  gradient 
flat  plate. 


Eckeimannv  IS 
Present  study 


V 


and  Eckelmann  (1970) 
r)  Intermi t tency  profil 


Figiirt*  .  Probability  density  function  of  the  longitudinal  turbulent 
velocity  P(u);  IJ  =9.36  m/sec,  R  =  3500. 


cycle  at  y*  -  3. A,  reported  by  Brodkey 


y,cm  y,cm 
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Figure  26. 


Comparison  of  mean-velocity  and 
profiles . 

a)  For  =  10.20  m/sec. 

b)  For  U*  =  8.55  ra/sec. 

00 


convective-veloci ty 


Figure  28a.  Peak  isocnrrelation  curves  from  the  correlation  of  the  boundary  layer  hot  wire 
and  the  surface  hot  wire  signals;  freestream  velocity  10.20  m/sec. 


als;  fri 


Repeatability  of  the  maximum  correlation  (R  ) 

^ T  u  max 
evaluation.  w 


ons  measured  at  several  lateral  positions  for 


z,cm  z,cm  2,cm 


Ratio  of  meon  surface  shear  with  and  without  the  plote  , 


height  above  the  surface. 


Probability,  P(t,,)  Probability,  P(t^) 


■’1 

1  ' 

'  1 
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Surface  shear,  N/m* 


Figure  37.  Probability  distribution  of  surface  shear  stress  due 
to  a  3  cm  wide  plate  placed  at  a  height  of  2  mm  above 
tlie  surface. 


138 


O 


Figure  38.  Probability  distributions  measured  with  a  thin  plate  placed  in  the  boundary 
layer  8  mm  (y*  ~  160)  above  the  surface. 


